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1.  Introduction 

A  large  body  of  the  literature  concerning  analytical  or  computa¬ 
tional  solutions  to  “impact”  problems  assumes  an  impact  condition 
in  the  form  of  an  applied  step  in  velocity  or  stress  on  the  face  of  the 
body  being  impacted  (Ma  and  Huang,  1996;  Lapczyk  et  al.,  1998; 
Wang  and  Sun,  2002;  Dayal  and  Bhattacharya,  2006;  Wuensche 
et  al.,  2009;  Sun  et  al„  2013;  Talebian  et  al„  2013;  Dutta  et  al„ 
2013).  However,  for  true  impact  in  which  at  least  one  of  the  bodies 
has  finite  thickness,  is  multilayered  or  is  continuously  inhomoge¬ 
neous,  multiple  reflected  or  scattered  waves  will  arrive  at  the 
impact  face  and  alter  the  stress  and  velocity  there.  Consequently, 
the  history  of  the  stress  and  the  particle  velocity  at  the  impact  face 
is  generally  time-dependent  and  may  be  difficult  to  determine  in 
advance  of  solution  to  the  problem. 

In  the  current  study,  we  consider  the  one-dimensional,  normal 
impact  problem  in  which  a  semi-infinite  flyer  collides  with  a  sta¬ 
tionary  target  plate  of  finite  thickness  i,  with  the  rear  face  of  the 
target  bonded  to  a  semi-infinite  backing  medium  (see  Fig.  1).  The 
flyer,  target  and  backing  medium  are  assumed  to  be  linear  elastic 
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and  homogeneous  with  generally  distinct  elastic  properties.  We 
also  assume  that  the  flyer  adheres  (or  “welds”)  to  the  front  face 
of  the  target,  in  the  sense  that  we  do  not  allow  separation  of  the 
flyer  and  target  after  impact  in  those  cases  where  wave  reflections 
would  otherwise  have  resulted  in  separation.1 *  Of  course,  this 
welding  assumption  is  superfluous  in  those  cases  where  separation 
would  not  have  occurred. 

We  are  primarily  interested  in  deriving  closed-form  expressions 
for  the  stress  and  particle  velocity  in  the  target  at  all  times  after 
impact.  The  solution  to  this  problem  is  simplified  by  applying 
appropriate  boundary  conditions  to  both  faces  of  the  target,  so  that 
the  one-dimensional  wave  equation  need  only  be  solved  in  the  tar¬ 
get  itself.  With  the  x-axis  oriented  as  shown  in  Fig.  1,  let  04  (x,  t) 
and  V\  (x,  t)  denote  the  stress  and  particle  velocity  at  the  point  x 
in  the  target  at  time  f,  and  let  z0  denote  the  impedance  of  the  flyer. 
Then  the  appropriate  impact  boundary  condition  on  the  front  face  of 
the  target  (x  =  0)  is 

z0Vi(0,t)±ai(0,t)=z0VoH(t),  (1) 


1  Such  a  physical  mechanism  is  known  to  govern  the  impact  restitution  of  colliding 

bodies  (Nicholas  and  Recht,  1990;  Shi,  1998;  Strange,  2000;  Goldsmith,  2001). 
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elastic  elastic  elastic 

flyer  target  half-space 


Fig.  1.  A  semi-infinite  flyer  impacting  a  target  oflength  /  with  impact  speed  V0.  The 
back  surface  of  the  target  is  bonded  to  a  half-space. 


where  H  denotes  the  Heaviside  step  function  and  V0  is  speed  of  the 
flyer  at  the  moment  of  impact  (t  =  0).  The  “+”  sign  in  (1)  is  used 
when  stress  is  taken  positive  in  compression  and  the  sign  when 
stress  is  taken  positive  in  tension.  Note  that  in  applying  (1),  the  par¬ 
ticle  velocity  and  stress  at  the  impact  face,  tq(0,  t)  and  or(0,  t),  are 
regarded  as  unknown— it  is  the  relationship  between  them  that 
determines  a  boundary  condition  corresponding  to  impact.  Further¬ 
more,  given  our  assumptions  on  the  flyer,  this  impact  boundary 
condition  is  (with  some  qualifications)  independent  of  the  configu¬ 
ration  and  properties  of  the  target  and  of  the  initial  conditions 
within  the  target  (see  Section  3.2).  Hence,  (1)  may  be  applied  to 
other  one-dimensional  impact  problems  where  interfacial  physics 
effects  such  as  adhesion  (Jayadeep  et  al.,  2014)  are  negligible. 
Alternatively,  for  computational  code  verification  purposes  (1)  is 
sufficient,  but  for  computational  code  validation,  one  may  also  need 
to  include  interface  physics  effects  (Jayadeep  et  al.,  2014; 
Israelachvili,  1991)  in  models  of  impact  phenomena. 

Condition  (1)  was  used  (without  proof)  by  Scheidler  and 
Gazonas  (2001)  in  a  study  of  impact  of  a  semi-infinite  flyer  on 
inhomogeneous  elastic  media;  see  also  Nicholas  and  Recht  (1990, 
Section  1.3.3)  for  a  different  impact  problem.  Chen  et  al.  (2004) 
have  emphasized  the  need  to  account  for  time  dependence  of  the 
stress  at  the  impact  face  for  plate  impact  on  periodically  layered 
media;  and  while  their  analysis  properly  accounted  for  this  time 
dependence,  they  did  not  state  or  use  the  impact  boundary  condi¬ 
tion  (1). 

Our  paper  is  organized  as  follows.  Section  2  contains  a  brief  dis¬ 
cussion  of  the  d’Alembert  solution  of  the  one-dimensional  wave 
equation  for  a  linear  elastic  solid;  this  is  limited  to  results  that  will 
be  used  in  the  sequel.  This  section  also  serves  to  establish  notation 
and  sign  conventions.  Section  3  contains  our  derivation  of  the 
impact  boundary  condition  (1 ),  which  utilizes  the  d’Alembert  solu¬ 
tion  of  the  wave  equation  in  the  flyer.  Some  alternative  forms  of 
this  condition  are  discussed,  and  we  determine  when  the  “impact” 
boundary  conditions  commonly  used  in  the  literature  are  correct  — 
either  exactly  or  in  some  asymptotic  sense. 

In  Section  4  we  derive  simple,  explicit  formulas  for  the  stress 
and  velocity  in  the  target  for  the  impact  problem  described  above. 
We  utilize  the  impact  boundary  condition  (1)  on  the  front  face  of 
the  target  and  an  appropriate  boundary  condition  on  the  back  face, 
so  that  the  one-dimensional  wave  equation  need  only  be  solved 
within  the  target.  We  show  that  the  stress  and  velocity  in  the  tar¬ 
get  can  be  obtained  from  the  solution  of  a  difference  equation  for 


one  of  the  d’Alembert  functions;  the  Laplace  transform  method  is 
used  to  obtain  this  solution  in  terms  of  a  series  of  Heaviside  func¬ 
tions.  Next,  from  the  Heaviside  series  we  obtain  simpler  expres¬ 
sions  for  the  two  d’Alembert  functions  in  terms  of  the  floor  (or 
greatest-integer)  function,  and  these  results  are  used  to  derive  sev¬ 
eral  simple,  explicit  formulas  for  the  stress  and  velocity  in  the  tar¬ 
get.  These  concise  relations  are  valid  for  all  times  after  impact  and 
all  points  in  the  target.  In  Section  5  we  use  these  results  to  evaluate 
the  velocity  and  stress  on  either  side  of  the  shock. 

Additional  properties  of  the  solutions  are  discussed  in  Sections 
6  and  7.  In  Section  6  we  determine  the  long  time  asymptotes  of  the 
velocity  and  stress  in  the  target  and  show  that  they  are  indepen¬ 
dent  of  the  properties  of  the  target.  We  also  determine  the  condi¬ 
tions  for  which  the  stress  at  the  impact  face  of  the  target  can 
(temporarily)  become  tensile  —  a  situation  that  indicates  separa¬ 
tion  of  the  flyer  and  target  would  have  occurred  were  it  not  for 
the  assumption  that  they  weld  on  impact.  In  Section  7  we  obtain 
solutions  for  the  special  cases  where  the  back  face  of  the  target 
is  either  stress-free  or  rigid.  These  cases  are  not  included  in  the 
main  results,  which  assume  that  the  impedance  of  the  backing 
medium,  z2,  is  nonzero  and  finite.  However,  solutions  for  the 
stress-free  and  rigid  cases  follow  easily  from  our  main  results  by 
taking  the  limits  as  z2  — >  0  or  oo,  respectively. 

Section  8  contains  plots  of  the  stress  and  velocity  histories  at 
either  the  front  face  or  the  midpoint  of  the  target  for  four  different 
impact  problems.  For  comparison  with  these  analytical  solutions, 
the  figures  also  include  the  discrete  solutions  of  a  recursive  system 
of  equations  for  impact  on  multi-layered  targets,  derived  else¬ 
where  (Gazonas  and  Velo,  in  preparation),  and  applied  to  the 
impact  boundary  value  problem  studied  here.  We  close  with  a 
discussion  of  various  methods  that  have  been  used  for  solving 
one-dimensional,  linear  elastic,  impact-type  problems  and  some 
conclusions  regarding  the  implications  of  the  techniques  used  here 
for  other  elastodynamic  problems  (Section  9).  Closed-form  expres¬ 
sions  for  the  jumps  in  stress  and  velocity  within  the  target  are 
derived  in  Appendix  A. 


2.  The  d’Alembert  solution  of  the  wave  equation 

This  section  contains  a  brief  discussion  of  the  d’Alembert  solu¬ 
tion  of  the  one-dimensional  wave  equation  for  a  linear  elastic  solid 
(Graff,  1975;  Eringen  and  §uhubi,  1975;  Achenbach,  1984; 
Drumheller,  1998;  Davison,  2008).  Only  those  results  used  in  the 
sequel  are  summarized  here.  This  section  also  serves  to  establish 
notation  and  sign  conventions.  The  subscript  k  can  take  the  values 
0,  1  or  2  for  the  flyer,  target  or  backing  medium,  respectively. 

2.1.  Solution  for  the  displacement 


We  consider  one-dimensional  deformations  of  a  homogeneous, 
linear  elastic  solid,  assuming  small  displacement  gradients.  Let 
uk(x,  t )  denote  the  displacement  of  material  k  along  the  direction 
of  impact  at  the  point  x  at  time  f,  measured  relative  to  the  initial 
configuration  at  t  =  0,  the  moment  of  impact.  Let  vk(x,  t)  and 
trfe(x,  t)  denote  the  corresponding  particle  velocity  and  the  stress: 


vk(x,  t) 


dUk 
dt  ’ 


ak(x,  t)  =  ±Ek 


dUk 
dx  ’ 


(2) 


where  Ek  is  the  elastic  modulus  and  duk/dx  is  the  longitudinal 
strain.  Here  and  below,  the  top  sign  corresponds  to  the  convention 
that  stress  is  positive  in  tension,  the  bottom  sign  to  the  convention 
that  stress  is  positive  in  compression. 

For  the  normal  plate  impact  problem  of  interest  here,  there  are 
no  other  displacement  or  strain  components.  The  analysis  in  this 
paper  can  also  be  applied  to  the  impact  of  isotropic,  linear  elastic 
rods  provided  that  lateral  inertia  can  be  neglected  (Graff,  1975); 
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the  appropriate  elastic  modulus  in  this  case  is  the  Young’s  modu¬ 
lus,  usually  denoted  by  E.  For  the  uniaxial  strain  state  in  normal 
plate  impact,  the  appropriate  elastic  modulus  is  the  longitudinal 
modulus.  For  isotropic  materials  this  is  given  by  X  +  2 p  or,  equiva¬ 
lently,  by  k  + 1  ft,  where  1,  ft  and  K  are  the  Lame,  shear  and  bulk 
moduli.  There  is  no  standard  notation  for  the  longitudinal  modu¬ 
lus,  so  we  have  simply  used  E  here. 

Momentum  balance  for  material  k  is  given  by 


,dok_  dvk 
±  dx  Pk  dt 


(3) 


where  pk  is  the  density,  assumed  positive.  The  relations  (2)  and  (3) 
yield  the  one-dimensional  wave  equation  for  the  displacement: 


2  cfuk  cfuk 
k  dx 2  “  dt2  ’ 


(4) 


where  ck  is  the  wave  speed.  Note  that  the  wave  equation  for  uk  is 
independent  of  the  sign  convention  for  the  stress.  d'Alembert’s 
general  solution  of  the  one-dimensional  wave  Eq.  (4)  is 


Mx,t)  =fk(t-x/ck)+gk(t  +  x/ck), 


(5) 


where  fk  andgfe  are  arbitrary,  sufficiently  smooth  functions  on  some 
interval  of  the  real  numbers  R.  Any  initial  or  boundary  conditions 
impose  restrictions  on  the  functions  fk  and  gk.  The  first  (second) 
term  on  the  right  in  (5)  represents  a  plane  longitudinal  wave  prop¬ 
agating  in  the  positive  (negative)  x-direction,  which  is  to  the  right 
(left)  with  the  materials  oriented  as  in  Fig.  1. 


2.3.  Conditions  at  the  shock  front 

In  the  remainder  of  the  paper  we  will  not  work  with  the  dis¬ 
placement  uk  directly.  Instead,  we  will  utilize  the  expressions  (6), 
(8)  and  (10)  for  the  velocity  vk,  stress  ok ,  and  d'Alembert  functions 
Fk  and  Gk  for  material  k.  For  the  impact  problem  considered  here, 
shock  waves  will  be  generated  in  the  target,  flyer  and  backing  med¬ 
ium.  Thus  uk  will  be  continuous  but  only  piecewise  smooth,  and 
vk,crk,Fk  and  Gk  will  be  piecewise  continuous  and  piecewise 
smooth  (in  fact,  piecewise  constant).  If  a  shock  front  is  located  at 
the  point  xs  in  material  k  at  time  ts,  then  vk(xs,ts)  and  ak(xs,ts) 
are  undefined,  at  least  physically,  whereas  their  one-sided  limits 
(see  Appendix  A)  exist  and  are  unequal,  by  definition.  On  the  other 
hand,  there  is  no  harm  in  assigning  a  value  to  the  velocity  and 
stress  at  the  shock  front  —  reasonable  choices  being  one  of  the 
one-sided  limits  or  their  average. 

Our  solutions  for  the  stress  and  velocity  in  the  target  will  be 
obtained  by  first  solving  for  the  d'Alembert  functions  Fj  and  G] 
and  then  using  (6)  and  (8)  to  obtain  V\  and  a i.  In  solving  for  the 
piecewise  continuous  functions  F i  and  Gi,  it  is  mathematically 
convenient  to  regard  these  functions  as  being  defined  at  all  points 
on  some  interval  of  R,  their  values  at  any  jump  discontinuities 
being  given  by  one  of  the  one-sided  limits.  Consequently,  our  solu¬ 
tions  for  V\  (x,  t)  and  <Ti  (x,  t)  will  have  values  assigned  on  the  shock 
front.  These  values  have  no  physical  significance,  and  the  use  of 
different  conventions  in  assigning  values  to  Ft  and  G i  at  their 
points  of  discontinuity  could  result  in  different  values  for  V\  and 
g \  on  the  shock  front. 


2.2.  Solution  for  the  velocity  and  stress 

Let  Fk  =  f'k  and  Gk  =  g'k,  where  a  prime  superscript  denotes  the 
derivative.  Then  by  (2)  and  (5),  the  velocity  and  stress  satisfy 

vk(x,  t)  =Fk(t-  x/ck)  +  Gk(t  +  x/ck),  (6) 

and 

£ 

ok{x,t)  =  ±-^[-Fk{t-x/ck)  +Gk(t  +  x/ck)].  (7) 

ck 

Henceforth,  we  take  stress  positive  in  compression,  so  that  the  bot¬ 
tom  (minus)  sign  holds  in  (2),  (3)  and  (7).  Then 


3.  The  impact  boundary  condition  for  a  semi-infinite  flyer 

We  begin  with  the  derivation  of  the  boundary  condition  on  the 
face  of  the  flyer,  using  the  d'Alembert  solution  of  the  wave  equa¬ 
tion  in  the  flyer.  This  result  may  also  be  derived  by  the  method 
of  characteristics.  The  corresponding  impact  boundary  condition 
on  the  target  is  obtained  in  Section  3.2.  An  alternative  form  of  this 
impact  boundary  condition  is  discussed  in  Section  3.3.  In  Sec¬ 
tion  3.4  we  determine  when  the  “impact”  boundary  conditions 
commonly  used  in  the  literature  are  correct  —  either  exactly  or 
in  some  asymptotic  sense. 

3.1.  Boundary  condition  on  the  flyer 


ak(x,  t)  =  zk[Fk(t  -  x/ck)  -  Gk(t  +  x/Cfc)],  (8) 

where  zk  is  the  impedance: 

£ 

Zk  =  PkCk  =~±=  y/pfc.  (9) 

Ck 

Unless  stated  otherwise,  we  assume  that  zk,pk,ck  and  Ek  are  positive 
and  finite  for  each  k  =  0, 1 , 2.  The  only  exceptions  occur  at  the  end 
of  Section  3.4  where  some  of  the  material  constants  in  the  flyer 
(k  =  0)  are  allowed  to  approach  zero  or  infinity,  and  in  Sections 
7.1  and  7.2  where  some  of  the  material  constants  in  the  backing 
medium  ( k  =  2)  are  allowed  to  approach  zero  or  infinity.  From  (6) 
and  (8),  we  see  that 

2 Fk(t  -  x/ck)  =  vk(x,  t)  +  ak(x,  t)/zk, 

2 Ck(t  +  x/ck)  =  vk(x,  t)  -  ak(x,  t)/zk.  '  J 

These  relations  are  useful  for  determining  the  restrictions  imposed 
on  Fk  and  Gk  by  initial  and  boundary  conditions.  The  functions  and 
(x,t)h->Fk(t-x/ck)  and  (x,  t)  i->  Gk(t  +  x/ck)  are  the  Riemann 
invariants;  they  are  constant  on  the  characteristic  curves 
t  -  x/ck  =  constant  and  t  +  x/ck  =  constant,  respectively. 


Consider  a  semi-infinite  flyer  which  is  initially  stress-free  and 
traveling  in  the  positive  x  direction  at  uniform  speed  V0  (see 
Fig.  1).  The  flyer  impacts  the  target  at  the  instant  t=0;  the 
target  and  the  half-space  backing  it  occupy  the  region  x  >  0,  and 
the  flyer  occupies  the  half-space  x  ^  0.  The  initial  conditions 
in  the  flyer  are 

v0(x,  0)  =  V0>0  and  <r0(x,  0)  =  0,  x  <  0.  (11) 

The  d’Alembert  solution  for  the  velocity  and  stress  in  the  flyer  is 
given  by  (6)  and  (8)  with  k  =  0.  On  setting  k  =  0  in  (10)i,  we  have 

Vo(x,t)  +  Go(x,t)/Zo  =  2Fo(t-x/c0),  x^O.tSsO.  (12) 

On  setting  t  =  0  in  (12)  and  using  the  initial  conditions  (11),  we  see 
that  V0  =  2 F0(-x/c0)  for  x  ^  0;  equivalently,  2 F0(t)  =  V0  for  t  >  0. 
Since  t  -  x/c0  >  0  for  x  and  t  as  in  (12),  it  follows  that 

v0(x,  t)  +  er0(x,  t)/z0  =  Vo,  x  <  0,  t  >  0.  (13) 

On  setting  x  =  0  in  the  above,  we  see  that  the  stress  and  velocity  at 
the  impact  face  of  the  flyer  are  related  by 

z0v0{0,t)  +  a0{0,t)=z0V0,  t  ^  0.  (14) 
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3.2.  Impact  boundary  condition  on  the  target 

Up  to  this  point  the  target  has  not  entered  the  analysis.  As  long 
as  the  flyer  and  target  are  in  contact,  their  displacements,  velocities 
and  stresses  must  be  continuous  across  the  impact  interface  and 
hence  coincide  at  the  interface.  Thus2 

ZA(0,t)  =  Vo(O,t),0i(O,t:)  =  <x0(0,t),  0<f<tc,  (15) 

where  a  “1”  subscript  is  used  for  target  variables  and  tc  is  the  dura¬ 
tion  of  contact  of  the  flyer  and  target,  which  may  be  finite  or  infi¬ 
nite.  We  do  not  assume  that  the  velocities  of  the  flyer  and  target 
coincide  at  f  =  0.  Indeed,  v0(0,0)  =  V0  by  (11),  whereas  the  target 
is  assumed  to  be  initially  at  rest  and  stress-free: 
zq(x,  0)  =  0,  (X,  (x,  0)  =  0,  0 Ux^l,  (16) 

where  1  is  the  length  of  the  target;  in  particular,  zq(0,0)  =  0. 
Thus  211(0,0)  ^  Vo(0.0),  and  more  significantly,  limxi0  V\  (x,  0)  ¥= 
limxtot,o(x,0). 

The  contact  relations  (15),  together  with  the  boundary 
condition  (14)  on  the  flyer,  yield  a  relation  between  the  stress 
and  particle  velocity  in  the  target  at  the  impact  face  x  =  0: 

ZoZ!i(0,t)  +  <7,(0,t)  =Z0V0,  0  <t^tc;  (17) 

see  also  Nicholas  and  Recht  (1990,  Eq.  (52))  for  impact  of  a  finite 
flyer  on  a  finite  target  (with  no  backing).  Our  derivation  of  (17) 
made  essential  use  of  the  one-dimensional  nature  of  the  impact 
problem  considered  here  and  also  of  the  assumptions  that  the  flyer 
is  homogeneous,3  linear  elastic  and  semi-infinite.  Since  no  informa¬ 
tion  about  the  target  was  utilized  in  the  derivation  of  (17),  it  follows 
that  this  impact  boundary  condition  is  essentially  independent  of  the 
configuration  and  properties  of  the  target  and  of  the  initial  conditions 
within  the  target.4  The  target  need  not  have  elastic  (let  alone  linear 
elastic)  response,  it  may  have  any  combination  of  discrete  or  contin¬ 
uously  varying  inhomogeneities,  and  it  may  be  initially  deformed 
and  undergoing  non-uniform  motion. 

If  the  duration  of  contact  is  finite,  then  the  front  face  of  the  tar¬ 
get  is  unstressed  once  separation  occurs,  and  the  appropriate 
boundary  condition  on  the  front  face  of  the  target  would  be 
ffi(0,  t)  =  0  for  tc  <  t  <  tr,  where  tr  is  the  first  instant  at  which 
the  flyer  and  target  regain  contact.  However,  the  duration  of  con¬ 
tact  tc  cannot  be  determined  in  advance  except  for  the  simplest 
problems.  These  complications  are  bypassed  by  assuming  (as  is 
done  here)  that  the  flyer  and  target  weld  on  impact.  Then  tc  =  oo, 
and  the  impact  boundary  condition  (17)  becomes5 

z0Vi{0,t)  +  <7i(0,  t)  =z0V0,  t>0.  (18) 


3.3.  Alternative  forms  of  the  impact  boundary  condition 

We  begin  by  considering  a  minor  variation  of  the  impact  bound¬ 
ary  condition  (18)  —  we  assume  that  it  also  holds  at  the  moment  of 
impact,  t  =  0: 


2  These  relations  do  not  strictly  hold  at  an  instant  t  =  ts  at  which  the  shock  front 
(reflected  from  the  back  surface  of  the  target)  arrives  at  the  flyer/target  interface, 
since  the  velocity  and  stress  are  undefined  at  this  instant.  However,  (15)  holds  in  the 
sense  that  zq(0,  t)  and  Vo(0,  ij  have  the  same  the  one-sided  limits  as  t  ts  from 
above  or  below,  and  similarly  for  the  stresses. 

3  It  can  be  shown  (by  a  different  method)  that  (17)  continues  to  hold  if  the  density, 
elastic  modulus  and  wave  speed  of  the  flyer  vary  with  position  in  such  a  way  that  the 
impedance  of  the  flyer  is  constant. 

4  This  statement  requires  some  qualifications.  The  duration  of  contact  tc  clearly 
depends  on  both  the  configuration  and  properties  of  the  target.  The  initial  velocity  of 
the  target  (if  any)  at  x  =  0  must  be  less  than  V0  in  order  for  impact  to  occur.  Any 
variations  in  target  properties  or  initial  conditions  must  be  in  the  x-direction  only,  so 
that  the  one-dimensional  nature  of  the  problem  is  retained  and,  in  particular,  so  that 
the  stress  and  particle  velocity  on  the  impact  face  are  uniform. 

5  If  stress  had  been  taken  positive  in  tension,  then  the  appropriate  boundary 
condition  would  be  zav i(0.  t)  -  <7i(0,  t)  =  z0V0. 


ZoZA(0,t)  +  t7,(0,t)  =  Z0V0,  f  >  0.  (19) 

Since  this  implies  Zot'i(O.O)  +  or(0, 0)  =  z0V0,  the  initial  conditions 
(16)  cannot  hold  at  x  =  0.  Instead,  the  initial  conditions  are 

V\ (x, 0)  =  0,  <7i(x,0)  =  0,  0  <x</,  (20) 


with  iq(0.0)  and  cri(O.O)  to  be  determined,  consistent  with  (19). 
The  impact  and  initial  conditions  (19)  and  (20)  differ  from  (18) 
and  (16)  only  at  a  single  point  in  (x,  t)  space,  namely  (x,t)  =  (0.0). 
One  would  not  expect  this  difference  to  appreciably  affect  the  solu¬ 
tion.  Indeed,  the  solutions  obtained  from  the  two  sets  of  impact  and 
initial  conditions  differ  only  in  the  values  assigned  to  the  stress  and 
velocity  on  the  shock  front;  cf.  the  discussion  in  Section  2.3. 

We  are  interested  in  determining  the  stress  and  velocity  in  the 
target  for  all  times  t  >  0.  Conditions  in  the  target  at  times  t  <  0  are 
irrelevant.  Nevertheless,  if  we  regard  V\  (x,  t)  and  or  (x,  t)  as  defined 
for  all  times  t  e  R,  then  we  are  free  to  assume  that  the  target  was 
stress-free  and  at  rest  for  all  times  prior  to  impact,  so  that 

Vi  (x.  t)  =  0.  or(x,  t)  =  0,  0  <  x  <  /,  t  <  0.  (21) 


Clearly,  the  boundary  condition  (19)  would  not  apply  for  t  <  0.  On 
the  other  hand,  if  H  denotes  the  Heaviside  step  function: 


(22) 


then  the  boundary  condition 

z0zq(0,t)  +  t7](0,  t)  =  z0V0H(t),  tsR  (23) 

is  consistent  with  both  the  impact  condition  (19)  and  with  (21). 
This  is  the  form  in  which  the  impact  boundary  condition  was  stated 
by  Scheidler  and  Gazonas  (2001).  The  conventions  described  in  this 
subsection  and,  in  particular,  the  impact  boundary  condition  in  the 
form  (23),  will  be  used  in  the  sequel. 

Finally,  recall  that  the  derivation  of  this  impact  boundary  condi¬ 
tion  makes  essential  use  of  the  assumptions  that  the  flyer  is  homo¬ 
geneous,  linear  elastic,  and  semi-infinite  and  that  the  flyer  adheres 
to  the  target.  Note  that  since  the  flyer  is  assumed  to  have  positive 
density,  it  necessarily  has  infinite  mass.  Consequently,  our  impact 
boundary  condition  is  not  appropriate  for  the  case  of  a  flyer  of 
finite  length  and  mass. 


3.4.  Common  "impact-type"  boundary  conditions 

The  two  most  common  "impact-type”  boundary  conditions 
used  in  the  literature  are  a  step  in  velocity  and  a  step  in  stress: 


zq(0,t)  = 

=  ViH(t), 

< 

V 

o 

t  6  R; 

(24) 

and 

ffi(O.t)  = 

2iH(t), 

M 

V 

o 

t  6  R. 

(25) 

Here  Vi  and  24  are  assigned  constants.  In  this  section  we  examine 
when  these  conditions  truly  describe  impact  —  either  exactly  or  in 
some  limiting  sense. 

First,  consider  the  case  where  both  the  flyer  and  the  target  are 
homogeneous,  linear  elastic  and  semi-infinite.  Then  both  (24)  and 
(25)  hold  with  Vi  and  Si  determined  by  the  impact  speed  V0  of 
the  flyer  and  the  impedances  z0  and  Zi  of  the  flyer  and  target: 

Vi=-^—V0,  Si  =ZiVi=^-V0:  (26) 

see  Davison  (2008,  Eq.  (3.21))  and  the  discussion  below.  Note  that 
Vi  <  V0.  It  is  easily  verified  that  the  relations  (24)-(26)  are  consis¬ 
tent  with  the  impact  boundary  condition  (23);  that  is,  if  (24)-(26) 
hold  then  so  does  (23).  Of  course,  when  the  assumptions  in  italics 
above  are  satisfied,  only  one  set  of  conditions,  either  (24)  and 
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(26)i  or  (25)  and  (26)2,  would  have  to  be  imposed  as  a  boundary 
condition  on  the  target.  From  (26)  it  follows  that 


Vo 

Zo 


—  zo+zi  V .  —  zo+zi  y, 

-  z,  -  £|  Z, 

V'o-V'i  1  ziV  o-£i  ’■ 


(27) 


Thus  for  a  given  flyer  impedance  z0,  the  boundary  condition  (24) 
holds  for  an  arbitrarily  prescribed  VG  provided  that  the  impact 
speed  of  the  flyer  is  given  by  the  first  expression  for  V0  in  (27); 
and  the  boundary  condition  (25)  holds  for  an  arbitrarily  prescribed 
Ei  provided  that  the  impact  speed  of  the  flyer  is  given  by  the  second 
expression  for  V0  in  (27).  Similarly,  for  a  given  flyer  speed  V0,  the 
boundary  condition  (24)  holds  for  an  arbitrarily  prescribed 
V]  <  V0  provided  that  the  impedance  of  the  flyer  is  given  by  the 
first  expression  forz0  in  (27);  and  the  boundary  condition  (25)  holds 
for  an  arbitrarily  prescribed  Ei  <  Zi  V0  provided  that  the  impedance 
of  the  flyer  is  given  by  the  second  expression  for  z0  in  (27). 

The  boundary  condition  (24)  or  (25)  is  often  used  for  “impact- 
type”  problems  in  which  the  conditions  in  italics  above  are  not  sat¬ 
isfied.  While  there  is  nothing  inherently  wrong  with  this,  neither 
(24)  nor  (25)  is  the  appropriate  boundary  condition  for  impact  in 
such  cases.  For  example,  consider  a  homogeneous,  linear  elastic 
target  of  finite  length  /,  with  or  without  a  backing  medium.  Let  t , 
denote  the  travel  time  in  the  target,  that  is,  the  time  it  takes  the 
shock  wave  to  travel  the  length  of  the  target: 

t.  =  '  •  (28) 

Cl 

Then  2 1,  is  the  time  it  takes  the  shock  wave  to  travel  from  the 
front  (impact)  face  of  the  target,  reflect  off  the  back  face,  and  return 
to  the  front  face.  If  (24)-(26)  are  indeed  valid  for  a  semi-infinite 
target,  then  one  would  expect  that  these  relations  hold  for  all  times 
t  <  2f»  in  the  finite  target,  but  generally  not  beyond  that.  In  Sec¬ 
tion  5.3  we  will  show  that  this  result  follows  easily  from  our  expli¬ 
cit  solutions.  For  a  semi-infinite  target,  the  relations  (24)-(26) 
follow  by  taking  the  limit  of  the  result  for  the  finite  target  as  its 
length  /,  and  hence  the  travel  time  t,,  approaches  infinity. 

Finally,  we  indicate  how  the  "impact-type"  boundary  conditions 
(24)  and  (25)  can  be  obtained  as  (different)  extreme  limiting  cases 
of  the  true  impact  boundary  condition  (23),  without  imposing  any 
conditions  on  the  target.  First,  let  V0  =  Vi  and  write  (23)  as 
V\(0.  t)  =  VjH(t)  -<7i(0 ,  t)/z0.  Then  (24)  holds  in  the  limit  as  the 
impedance  z0  of  the  flyer  becomes  infinite.  By  (9)  with  k  =  0,  this 
latter  condition  corresponds  to  a  flyer  of  infinite  modulus  and  non¬ 
zero  density  (i.e.,  a  rigid  flyer),  or  a  flyer  of  infinite  density  and  non¬ 
zero  modulus.  Next,  let  V0  =  Ei/z0,  in  which  case  (23)  can  be 
written  as  a^(O.t)  =  Z,H(t)  -  z0Vi(0.t).  Then  we  obtain  (25)  in 
the  limit  as  the  impedance  of  the  flyer  goes  to  zero  with  Ei  fixed, 
in  which  case  the  speed  of  the  flyer  becomes  infinite.  By  Eq.  (9),  the 
flyer  has  zero  impedance  if  the  density  is  zero  and  the  elastic  mod¬ 
ulus  is  positive  and  finite,  or  vice  versa. 


4.  Solution  in  the  target  by  combined  d’Alembert  and  Laplace 
transform  methods 

In  this  section  we  derive  explicit  formulas  for  the  stress  and 
velocity  in  the  target  for  the  impact  problem  described  in  the  intro¬ 
duction.  We  utilize  the  impact  boundary  condition  on  the  front 
face  of  the  target  (Section  3)  and  an  appropriate  boundary  condi¬ 
tion  on  the  back  face  (Section  4.1).  In  Section  4.2  we  express  the 
impact  problem  in  non-dimensional  form  and  show  that  the 
non-dimensional  stress  and  velocity  in  the  target  can  be  obtained 
from  the  solution  of  a  difference  equation  for  one  of  the  non- 
dimensional  d’Alembert  functions.  In  Section  4.3  we  use  the 
Laplace  transform  method  to  obtain  the  solution  of  this  difference 
equation  in  terms  of  a  series  of  Heaviside  functions.  A  simpler  form 


of  the  solution  for  the  two  non-dimensional  d'Alembert  functions 
in  terms  of  the  floor  (or  greatest-integer)  function  is  obtained  in 
Section  4.4,  and  these  results  are  used  in  Section  4.5  to  obtain  sev¬ 
eral  simple,  explicit  formulas  for  the  non-dimensional  stress  and 
velocity  in  the  target.  These  results  are  recast  in  dimensional  form 
in  Section  4.6.  We  follow  the  conventions  in  Section  3.3;  in 
particular,  we  regard  the  stress  and  velocity  as  defined  for  all  times 
t  e  R.  However,  the  explicit  formulas  for  the  stress  and  velocity  in 
Section  4.6  are  only  valid  for  t  ^  0. 

4.1.  Back  face  boundary  condition 

The  d’Alembert  solution  for  the  velocity  and  stress  in  the 
half-space  backing  the  target  are  given  by  (6)  and  (8)  with  k  =  2. 
On  setting  k  =  2  in  (10)2,  we  have 

2G2(t  +  x/c2)  =  v2(x,t)  -  a2(x,t)/Z2,  t  e  R,  x  > /.  (29) 

The  domain  of  G2  is  all  of  R.  Since  the  half-space  is  at  rest  and 
unstressed  for  all  times  t  0,  the  right-hand  side  of  (29)  is  zero 
for  t  sg  0  and  x  >  /.  As  x  and  f  range  over  these  semi-infinite 
intervals,  t  +  x/c2  takes  on  all  real  values.  Hence  C2  =  0,  and  (29) 
reduces  to 

<r2(x,  t)  =  z2v2(x,  t),  t  e  R,  x  >  /.  (30) 

In  particular,  at  the  interface  with  the  target  (x  =  /),  the  velocity  and 
stress  in  the  half-space  are  related  by 

a2(l,t)  =z2v2(l,t),  teR.  (31) 

Then  by  continuity  of  stress  and  velocity  across  material  interfaces, 
the  velocity  and  stress  in  the  target  at  x  =  /  are  related  by6 

G\(l,t)  =z2Vi(l,t),  teR.  (32) 

This  boundary  condition,  like  the  impact  boundary  condition  (18),  is 
independent  of  the  properties  of  the  target. 


4.2.  The  d’Alembert  solution  in  the  target 


From  the  d’Alembert  solution  of  the  wave  equation,  the  velocity 
and  stress  in  the  target  are  given  by  (6)  and  (8)  with  k  =  1: 

V\  (x,  t)  =  Fj(t-  x/c, )  +  Gi  (t  +  x/c, ) , 

O’,  (x,  t)  =z1[F1(t-x/c1)  -  Gi  (t  +  x/c, )] . 


for  t  e  R  and  0  ^  x  ^  l.  The  domains  of  F,  and  G i  are  all  of  R.  Recall 
that  the  initial  conditions  are  given  by  (20)  and  the  conditions  prior 
to  impact  by  (21).  On  setting  k  =  1  in  (10)  and  using  (20)  and  (21), 
we  obtain  the  following  restrictions  on  the  d’Alembert  functions  F , 
and  Gi : 

F, (t)  =  0,  t  <  0  and  G,(t)  =  0,  t  <  t,  =  l/cx.  (34) 


The  derivation  of  the  solutions  for  F ,  and  G i,  and  the  subsequent 
expressions  for  V\  (x,  t)  and  <Ti(x,  t),  are  simplified  if  we  cast 
the  problem  in  non-dimensional  form.  We  define  the 
non-dimensional  distance  y  and  time  t  by 


Clt 

l  ’ 


(35) 


then  0  ^  y  ^  1  and  tsB.  For  t  >  0,  Ci  t  is  the  total  distance  (includ¬ 
ing  reflections)  that  the  wave  has  traveled  from  the  moment  of 
impact;  hence  z  is  also  the  total  non-dimensional  distance  the  wave 
has  traveled  from  the  moment  of  impact.  The  velocity  and  stress  are 
often  non-dimensionalized  by  dividing  by  the  wave  speed  and  elas- 


6  If  stress  had  been  taken  positive  in  tension,  then  the  appropriate  boundary 
condition  would  be  r 7,(1, f)  +  z2zq(/,t)  =  0;  see  also  Nonaka  et  at  (1996,  Eq.  (17)). 
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tic  modulus,  respectively.  For  the  impact  problem  considered  here,  Also  note  that  by  (47)  and  (46)i,  the  condition  (46)2  can  be 
it  turns  out  to  be  more  convenient  to  use  strengthened  to 


v{y,  x)  = 


V\  (x,t) 
V,  ' 


a(y,  x)  = 


<Ti(x,t) 


(36) 


for  the  non-dimensional  velocity  and  stress.  Recall  that  VG  and  Si 
are  given  in  terms  of  the  impact  speed  V0  and  the  impedances  z0 
and  Z\  of  the  flyer  and  target  by  the  expressions  in  (26).  This  choice 
is  motivated  by  the  fact  that  V i  and  Si  are  the  initial  jumps  in 
velocity  and  stress  in  the  target  (see  Section  5);  however,  at  this 
point  we  simply  regard  the  expressions  in  (26)  as  definitions  of 
V\  and  Si. 

On  dividing  the  impact  boundary  condition  (23)  by  z0V0  and 
then  using  (35),  (36)  with  x  =  0,  the  expressions  (26)  for  V\  and 
Si,  and  the  fact  that  H(t)  =  H( r),  we  obtain 


z0 

Zp+Zi 


v(0,x) ■ 


Zo+Zi 


a(0,x)  =  H(x),  xe 


(37) 


Similarly,  on  dividing  the  back  face  boundary  condition  (32)  by  Si 
and  using  (26),  (35)  and  (36)  with  x  =  /,  we  obtain 


<T(l,x)  =  fv(i,x), 


t  e 


(38) 


These  are  the  non-dimensional  forms  of  the  boundary  conditions  on 
the  target.  Since  the  d’Alembert  functions  F i  and  Gi  have  the 
dimensions  of  velocity  and  since,  by  (35), 


t±x/c,  =  r.(T±y), 


(39) 


we  define  non-dimensional  d’Alembert  functions  F  and  G  as 
follows: 


F(  x)  = 


Fi(t.x) 

V, 


C(x)  = 


G,(t. r) 

V, 


(40) 


Then  by  (33),  (36),  (39)  and  (40),  the  d’Alembert  solution  for  the 
non-dimensional  velocity  and  stress  in  the  target  is 


v(y,x)  =F(x-y)  +  G(T+y), 
<7(y,T)=F(r-y)-G(T  +  y), 


(41) 


for  t  6  R  and  Osjy^  1.  The  domains  of  F  and  G  are  all  of  R. 

The  following  non-dimensional  parameters  occur  repeatedly  in 
the  sequel: 


a  = 


Zp-Zj 
Zp  +  zC 


z2  -  z  i 
'  z2  +  zr 


A  =  ab. 


(42) 


G(x)  =  0,  x  <  2.  (49) 

In  the  next  section,  we  will  see  that  the  difference  Eq.  (48)  for  F, 
subject  to  the  initial  condition  (46)i,  is  easily  solved  by  the  Laplace 
transform  method.  Then  the  solution  for  G  follows  from  (47),  and 
the  non-dimensional  velocity  and  stress  can  be  obtained  from  (41 ). 

4.3.  Solution  for  F  by  Laplace  transforms 


Let  /  denote  a  real-valued,  piecewise-continuous  function 
defined  on  R,  with 

/( T)  =  0,  x  <  0.  (50) 

We  employ  the  right-sided  or  unilateral  Laplace  transform  of  / 

(Doetsch,  1974): 

poo 

/(s)  =  G{/}(s)=Gt{/(t)}(s)=  /  f(x) e~STdx,  (51) 

Jo 

where  s  =  a  +  ioj  is  generally  a  complex  number  with  91(s)  >  0.  The 
Laplace  transform  of  the  Heaviside  step  function  is  given  by 

H(s)=I.  (52) 

We  also  need  the  Laplace  shifting  theorem  for  functions  satisfying 

(50):  if  t,  >  0,  then 

Cx{f(x-Xi)}(s)  =  e-^f(s).  (53) 

On  taking  the  Laplace  transform  of  the  difference  Eq.  (48)  for  F  and 
using  (52),  (46)!,  and  (53)  with /  =  F  and  Ti  =  2,  we  obtain 

(l-Ae-2s)F(s)  =  l.  (54) 


Since  9?(s)  >  0,  we  have  |e_2s|  =  e~29i(s)  <  1,  so  that  by  (43  )3, 
[Ae~2s|  <  |A|  <  1.  Hence  1  -  Ae~2s  ^  0,  and  (54)  yields  the  following 
expressions  for  the  Laplace  transform  of  F: 


F(s)  = 


1 


1 


s  1  -Ae~ 


l+£(Ae- 


-i  oo 


.1 


(55) 


The  inverse  transform  of  (55)  yields  the  non-dimensional  d'Alem¬ 
bert  function  F.  On  using  (52)  and  the  shifting  theorem  (53)  with 
/  =  H  and  Ti  =  2n,  we  see  that 


Since  the  impedances  are  positive, 

|a|<l.  \b\  <  1,  |A|  <  1 .  (43) 

On  substituting  (41)  with  y  =  0  into  the  non-dimensional  impact 
boundary  condition  (37)  and  rearranging  terms,  we  obtain 

F(t)  +  aG(x)  =  H(x),  xeR.  (44) 


1  e~2ns  =  jCz{H(x  -  2n)}(s).  (56) 

Then  by  (56)  and  (52)  it  follows  that  F  is  given  by 

oo 

F(x)  =  H(x)  +  Y)AnH(x-2n),  xeR;  (57) 

n= 1 


Similarly,  on  substituting  (41 )  with  y  =  1  into  the  non-dimensional 
back  face  boundary  condition  (38)  and  rearranging  terms,  we  obtain 

G(t  +  1)  =  —bF(x  —  1),  T  e  R.  (45) 

Also,  by  (40),  the  non-dimensional  form  of  the  initial  conditions 

(34)  is 

F(t)  =  0,  t  <  0  and  G(t)  =  0,  t<1.  (46) 

Next,  observe  that  the  relation  (45)  is  equivalent  to 

G(x)  =  -bF(  r  —  2),  xeR.  (47) 

On  substituting  this  relation  into  (44)  and  using  (42  )3,  we  obtain  the 
following  difference  equation  for  F: 

F(t)  -  AF(t  -  2)  =  H(x),  xeR.  (48) 


and  by  (47), 

oo 

G(x)  =  —bH(x  -  2)  —  b^A"H(T  -  2(n  +  1)),  xeR.  (58) 

n-l 

It  follows  that  F  has  jump  discontinuities  at  the  non-dimensional 
times  r  =  0,2,4,...;  and  G  has  jump  discontinuities  at  the  non- 

dimensional  times  t  =  2.4,6 _ 

If  x  <  0,  then  H(x)  =  FL(x  -  2 n)  =  0,  and  we  recover  the  initial 
condition  (46)3  from  (57).  More  generally,  for  any  tell  there  are 
only  a  finite  number  (possibly  zero)  of  nonzero  terms  in  the  sum 
in  (57),  since 

H( x  -  2 n)  =  0  iff  x  -  2n  <  0  iff  n  >  t/2;  (59) 

and  FI(x-2n)  =  l  otherwise,  that  is  for  n^x/2.  Thus  (57) 
reduces  to 
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F(t)  =  H(x)  +  J2  A",  T  e  OS,  (60) 

l^n^t/2 

provided  the  sum  is  regarded  as  zero  if  there  is  no  integer  n  satisfy¬ 
ing  1  <  n  <  r/2,  that  is,  if  r/2  <  1.  An  analogous  result  holds  for  C. 

4.4.  Solution  for  F  and  G  in  terms  of  the  floor  function 


The  concise  expression  in  (68)  is  not  valid  for  r  <  0  since  it  yields  a 
nonzero  result,  contrary  to  (49).  However,  for  0  <  t  <  2  we  have 
\x/2\  =  0  and  thus  ALt/2J  =  1,  so  that  G(x)  =  0,  consistent  with 
(49).  For  2  sc  t  <  4  we  have  [t/2J  =  1,  so  that  (68)  yields 

G(r)  =  -fa,  2  <  t  <  4.  (69) 


If  r  is  not  a  positive  even  integer,  then  there  will  be  no  integer 
n  ^  1  that  is  equal  to  r/2,  in  which  case  it  would  not  be  correct  to 
write  the  sum  in  (60)  as  Let  Z  denote  the  set  of  integers. 

The  floor  function  (or  greatest-integer  function)  |_-J  :  R  — >  Z  is 
defined  by  the  condition  that  [rj  is  the  greatest  integer  less  than 
or  equal  to  the  real  number  r.  Some  properties  of  the  floor  function 
are  listed  below;  here  r  is  a  real  number  and  n  and  m  are  integers: 

n  <  r  <  n  +  1  LrJ  =  n-  Lr±nJ  =  LrJ±n>  (61) 

n  <  r  «=>  n  ^  LrJ  i  r  <m  |_rJ  ^  171  ~  1  ■  (62) 


4.5.  Expressions  for  the  non-dimensional  velocity  and  stress  in  the 
target 

On  substituting  the  expressions  (65)  and  (68)  for  F  and  G  into 
the  relations  (41)  for  the  non-dimensional  velocity  and  stress,  we 
obtain  the  compact  formulas 

v(y, x)  =  y-Lj  [l  -  AL(T_y)/2J+1  -  fa(l  -  A^+y>/2J)] ,  (70) 

ff(y,T)  [l  _Al(T~y)/2J+1  +b(l  _AL(T+y)/2J)].  (71) 


On  setting  r  =  t/2  in  the  left  relation  in  (62),  we  see  that  (60)  can  be 
written  as 


b/2J 

F{x)  =  H(x)  +  J2A"  ’  TeR’ 

n=l 


(63) 


provided  the  sum  is  regarded  as  zero  if  [t/2J  <  1,  equivalently,  if 
t  <  2.  In  particular,  (63)  yields 

{0,  T  <  0, 

1.  0  s;  T  <  2,  (64) 

1+A,  2s;t<4. 

Next,  we  claim  that 
i  _  /ib/2j+i 

F(t)=— j— ,  t>-2,  (65) 


provided  that  we  use  the  convention  0°  =  1  when  A  =  0.  Since 

A0  =  1  for  any  nonzero  A,  this  convention  implies 

A°=l.  (66) 

for  all  A.  To  prove  the  claim,  first  note  that  for  any  real  number  q^l 
and  any  integer  N  >  1, 


N 


!+E  9' 

n=l 


1  -  qN+1 
1  -q 


(67) 


This  relation  also  holds  for  N  =  0  provided  the  sum  is  regarded  as 
zero  in  this  case.  Suppose  that  x  >  0.  Then  by  (63), 
F(t)  =  1  +  5Zj,I/12JAn.  Since  A^l  and  since  \x/2\  0  in  this  case, 
we  may  apply  (67)  with  q=A  and  N  =  [t/2J,  which  yields  the 
expression  in  (65).  Now  suppose  that  -2  sg  x  <  0.  In  this  case 
|t/2J=-1,  so  the  right-hand  side  of  (65)  reduces  to 
(1  -  A°)/(l  -  A),  which  is  zero  by  the  convention  (66).  Since 
F(x)  =  0  for  x  <  0,  (65)  gives  the  correct  result  for  -2  sg  x  <  0, 
which  proves  the  claim.The  concise  expression  in  (65)  is  not  valid 
for  t  <  -2  since  it  yields  a  nonzero  result.  However,  a  closed- 
form  expression  for  F(r)  for  t  >  -2  is  all  that  is  really  needed  here, 
since  by  (47)  it  yields  a  closed-form  expression  for  G(t)  for  x  >  0. 
Indeed,  since  |_(t-2)/2J  =  (_t/2  —  1J  =  |r/2J  -  1,  from  (47)  and 
(65)  we  obtain 

1  -  A[r/2J 

G(x)  =  -b  1  _  A  ■  ?>0.  (68) 


These  relations  and  those  below  are  valid  for  all  0  sg  y  sg  1  but 
only  for  x  >  0.  At  a  given  non-dimensional  position  y  in  the 
target,  v(y ,  t)  and  a(y,  t)  suffer  jump  discontinuities  at  those 
non-dimensional  instants  x  ^  0  for  which  either  L(t  -  y) /2J  or 
L(r  +  y)/2J  jump  to  the  next  integer  value. 

To  shorten  subsequent  expressions,  we  will  use  an  abbreviated 
notation  for  the  exponents  in  (70)  and  (71).  By  (28)  and  (35)  (see 
also  (39)),  we  have 
T±y_t±x/Ci  _c^t±x 

2  "  2 1,  21  '  (  ' 


Let 


M(y,T)  = 


it+yi 

t  +  x/c, 

Cit  +  X 

L  2  J 

21, 

21 

=  M(x,  t), 


(73) 


and 

N(y,T)=[^J+l 


t  -  x/c, 
2 1, 


+  1 


Cit  -  X 
~ 21- 


+  1  =  N(x,  t). 


(74) 


M(y,  r)  and  N(y,  t)  are  piecewise-constant,  integer-valued  functions 
ofy  and  x  which,  for  fixed  y,  are  nondecreasing  functions  of  x.  While 
these  functions  are  well-defined  for  all  real  y  and  t,  for  use  in  (70) 
and  (71)  we  consider  only  their  values  for  x  >  0  and  0  <y<i- 
On  rearranging  the  terms  in  (70)  and  (71)  and  using  (73)  and 
(74),  we  obtain 

v (y, T)  =  j-E  ^l  -b- AN(y,T)  +  faAM(y,T)] ,  (75) 


a (y,  t)  =  y-E  [l  +  fa  -  AN(y'T)  -  faAM(y’T)] . 

From  (42)  we  have 

^  _  zo  ~  Z1  z2  ~  zi  ^  _  z2  ~  Z1 
~Z0+Z,  Zj+Zj  ’  +Z]  ' 

It  follows  that 

1  ^  m  2zi  z0  +  z2 

_  Zi  +  Z2  Zo  +  z,  ’ 

i-fa  =  EZi_  1  +  b  =  EZ2_. 

Zi  +  z2  Z]  +  z2 


(76) 


(77) 


(78) 


4.6.  Expressions  for  the  velocity  and  stress  in  the  target 


7  The  ceiling  function  [•]  :  R  — »  Z  is  defined  by  the  condition  that  \r]  is  the  smallest 
integer  greater  than  or  equal  to  the  real  number  r.  Any  expression  involving  floor 
functions,  such  as  (63),  can  be  converted  to  an  equivalent  expression  in  terms  of 
ceiling  functions  by  using  the  floor-to-ceiling  transformation  [rj  =  —  r]. 


From  the  non-dimensionalization  (36),  we  see  that  the  particle 
velocity  V\  and  stress  are  related  to  the  non-dimensional  veloc¬ 
ity  and  stress,  v  and  a,  by 
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Vi(x,t)  =  Viv(y,  x),  <Ti  (x,  t)  =  X,  a  (y,  r) . 


(79)  5.1.  Evaluation  of  the  exponents  M(x ,  t)  and  N(x,  t) 


Expressions  for  zq(x,  t)  and  a^(x,t)  for  0  <  x  sg  /  and  t  3=  0  follow 
from  (79),  the  relations  for  v(y,x)  and  o(y,  x)  in  Section  4.5,  and 
the  fact  that  M(y,  r)  =  M(x,  t)  and  N(y,  z)  =  N(x,  t)  (see  (73)  and 
(74)).  For  example,  from  (75)  and  (76)  we  obtain 


vl(x,t)  =  ^[\-b-AN™  +  bAMM], 


a,  (x.  t)  =  ^  [l  +  b  -  Amt)  -  Mmw)]  . 


(80) 


(81) 


The  piecewise-constant,  integer-valued  functions  M(x,  t)  and  N(x,t) 
are  nondecreasing  functions  of  t  for  fixed  x.  Additional  properties  of 
these  functions  are  derived  in  Section  5.1.  For  comparison  with  the 
concise  expressions  (80),  (81),  we  state  the  solutions  in  the  more 
traditional  form  of  an  infinite  series  of  Heaviside  functions: 


vi(x,  t)  =  V1Y)An[H{t  -  tj)  -  bH(t  -  t2)], 


ff,(x,t)  =  X,  ^A"[H(t  -  t,)  +  bH(t  -  t2)]. 


(82) 


(83) 


ti  =- 


(84) 


where 

2ln  +  x  f  2/(n+l)-x 

Cl  C] 

These  can  be  obtained  by  substituting  the  series  solutions  (57)  and 
(58)  for  F  and  G  into  (41)  and  using  (79). 

Recall  that  A  and  b  are  given  by  (77),  1  ±  b  and  1  -  A  by  (78), 
and  Vi  and  X!  by  (26).  Alternative  expressions  for  the  leading  coef¬ 
ficients  in  (80)  and  (81)  are  given  by 

Vi  z0z,  +z2V0 


1  -  A  Zi  z0  +  z2  2 


2i 


V, 


Zq 


Zi  +  z2  v0 


(85) 


1-A  1-A  z0+z22' 

For  the  next  set  of  results,  it  is  convenient  to  introduce  the  terms 


V  = 

v  oo  — 


z  0 


Zo+Z2 


i/0=!0+!'  v,. 


—  Z2 17^  — 


Zo+Z2 


Z0Z2  ..  z0+z,  z2 

-  Vo  = - 2-] , 

z0  +  z2  z0  +  z2  z, 


(86) 


(87) 


where  the  expressions  on  the  right  follow  from  (26).  Clearly,  and 
Xoo  are  positive:  in  Section  6.1  we  will  show  that  they  are  the  long 
time  asymptotes  of  velocity  and  stress  in  the  target.  On  substitution 

of  (77)2,  (78)2.3,  (85)— (87)  into  (80)  and  (81),  we  obtain  the  concise 
expressions 


V\{x,  t)  =  ^p 


(Ti(X,t)  =- 


:  -  1)  A"M  +  (r~1)  aM<V> 


Z2 


2-  ( -+i)an(v)  +  (--  1 


z2 


(88) 


(89) 


We  begin  by  determining  the  subsets  of  the  (x,  t)  plane  on 
which  the  integer-valued  functions  M(x,  t)  and  N(x,  t)  are  constant. 
Let  m  be  any  integer.  From  the  definition  (73),  we  see  that 


M(x.  I)  =  m 


t  +  X/Ci 


20 


=  m 


t  +  X/Ci 

m^  '  <  m  +  1 


2mt»  -  x/ci  <  t  <  2 (m  +  l)t»  -  x/c  1 

Ci(2mt,  -  t)  <  x  <  Ci[2(m  +  1)0  -  t]. 
From  the  definition  (74),  we  see  that 
N(x,  t)  =  m 


t  -  X/C! 


20 


=  m  -  1 


.  t  -  X/Ci 

<;=>  m-  1  ^  ^ —  <  m 


2{m  -  1)0  +  x/ci  ^  t  <  2m0  +x/Ci 
Ci(t  -  2m0)  <  x  <  Ci[t  —  2(m  -  1)0]- 


(90a) 

(90b) 

(90c) 

(90d) 

(90e) 

(91a) 

(91b) 

(91c) 

(91  d) 
(91e) 


To  simplify  the  subsequent  discussion  of  these  results,  we  introduce 
abbreviations  for  some  of  the  expressions  above  (the  arrow  super¬ 
scripts  will  be  motivated  in  the  next  subsection): 


s^(x)  =  2m0  -  x/ci  =  (2m  -  x//)0, 

s^(x)  =  2 (m  -  l)0+x/Ci  =  [2(m  -  1)  +x/(]0 

and 

S;(t)  ee  Ci(2m0  -  f)  =  2ml- c^t, 

S^(t)  =  Ci  [t  -  2 (m  -  1)0]  =  c,  t  -  2(m  -  1)1. 


Observe  that 

t  =  4W 
t  =  4W 


x  =  S„(t), 

■  X  =  Sm( f)- 


(92) 


(93) 


(94) 


Using  (92)  and  (93),  we  can  rewrite  the  bottom  two  equivalences  in 

(90)  and  (91)  as 


M(x,  t)  =  m 
N(x,  t)  =  m 


'  Sm  (X)  ^  t<4+l(X), 

•S;(t)^x<S;+1(t), 

■Sm(X) 

■  Sm+A[)  <  X  <  Sm(t)- 


(95) 


5.  Solutions  on  either  side  of  the  shock 


The  relations  (90)-(95)  are  valid  for  any  real  numbers  x  and  t.  How¬ 
ever,  the  expressions  (80),  (81)  and  (88),  (89)  for  the  velocity  and 
stress  in  the  target  are  only  valid  for 


The  concise  relations  for  the  velocity  and  stress  in  Section  4.6 
are  valid  for  all  times  t  >  0  and  all  points  0  ^  x  ^  1  in  the  target. 
We  begin  this  section  by  evaluating  the  exponents  lW(x,  t)  and 
N(x,  t)  in  these  relations.  These  results  are  used  in  5.2  to  determine 
the  position  of  the  shock  at  each  instant  and  in  Section  5.3  to  eval¬ 
uate  the  velocity  and  stress  on  either  side  of  the  shock.  In  Appendix 
A  we  use  the  results  in  Sections  3-5  to  derive  several  expressions 
for  the  jumps  in  velocity  and  stress  across  the  shock. 


0  <  x  <  1  and  t  >  0.  (96) 

These  restrictions  are  assumed  in  the  sequel.  A  graphical  represen¬ 
tation  of  the  inequalities  (95)  and  (96)  is  given  in  Figs.  2  and  3. 

Subject  to  the  restrictions  (96),  for  m  =  1,2,3,...  the  set  of 
points  in  the  (x,  t)  plane  given  by  t  =  s^(x)  is  the  mth  line  of  jump 
discontinuity  of  the  function  M;  similarly,  t  =  s^(x)  is  the  mth  line 
of  jump  discontinuity  of  the  function  N.  These  lines  intersect  only 
at  x  =  0  and  x  =  I,  that  is,  at  the  front  and  back  faces  of  the  target. 
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For  x  =  0  this  occurs  at  t  =  2t,.  4t„,  6t», . . for  x  =  /  the  intersec¬ 
tions  occur  at  t  =  t, ,  3t, ,  5t„ , . . .. 

5.2.  The  shock  wave 

Since  |Aj  <  1,  it  follows  that  and  AN(x,t>  suffer  jump  discon¬ 
tinuities  along  the  same  lines  in  the  (x,  t )  plane  as  the  exponents  M 
and  N.  From  the  relations  (80)  and  (81)  or  (88)  and  (89),  these  are 
also  the  lines  of  jump  discontinuity  of  V\  (x,t)  and  <7i(x,  t)  in  the 
(x,  t)  plane.  In  other  words,  the  locus  of  the  shock  in  the  (x,  t)  plane 
is  the  union  of  the  lines  t  =  s^fx)  and  t  =  s^(x)  (equivalently,  the 
lines  x  =  S^(t)  and  x  =  S^(t)),  subject  to  the  restrictions  (96),  as 
m  ranges  over  all  positive  integers.  This  is  illustrated  in  Fig.  4.8 9 

The  line  x  =  S^(t)  has  positive  slope  and  represents  the  mth 
forward  propagating  shock  (m  =  1,2,3,...).  Note  that  the  forward 
direction  is  to  the  right  when  the  target  is  oriented  as  in  Fig.  1 ; 
hence  the  right  arrows.  The  line  x  =  S^(t)  has  negative  slope  and 

represents  the  mth  backward  propagating  shock  (m  =  1,2,3 _ ). 

The  backward  direction  is  to  the  left  when  the  target  is 
oriented  as  in  Fig.  1;  hence  the  left  arrows.  The  changes  in  the 
direction  of  propagation  of  the  shock  at  t  =  t,,3t„,  5f„, . . .  and 
f  =  2t„,4f,,6f„,...  are  due  to  reflection  from  the  back  and  front 
faces  of  the  target,  respectively.  The  triangular  shaped  regions  in 
Fig.  4  are  the  regions  in  which  both  M  and  N  are  constant.  The 
relations 


Flyer  f  Target 

2(m+l)t * 

(2m+l)tt 

2m  t » 

(2m-l)tt 
2(m-l)t *  — 


t*  -  ■ 

M(x,t)  =  0 

0 - - 

0  1  x 

Fig.  2.  Regions  of  the  (x,t)  plane  in  which  the  integer-valued  functions  M  are 
constant.  The  values  indicated  are  assumed  on  the  solid  lines  but  not  on  the  dashed 
lines. 


N(x,t)>M(x,t)>0, 

N(x,t)=M(x,t )  or  M(x,  t)  +  1. 

are  evident  in  this  figure  and  also  in  Figs.  2  and  3;  analogous  rela¬ 
tions  hold  for  the  functions  M(y,T)  and  N(y,  r)  in  (73)  and  (74).  By 
(80)  and  (81)  or  (88)  and  (89),  it  follows  that  the  triangular  regions 
in  Fig.  4  are  also  regions  of  constant  velocity  and  stress.  As  indicated 
in  the  figure,  V*,  VR.  VR, . . .  and  £jj,  E*,  E*, . . .  denote  the  sequences 
of  (constant)  velocity  and  stress  values  to  the  right  of  the  shock; 
V\ .V'2.  V'3.  . . .  and  E^EjjE?;,...  denote  the  sequences  of  (constant) 
velocity  and  stress  values  to  the  left  of  the  shock.  Simple  expres¬ 
sions  for  these  values  will  be  derived  in  the  next  subsection.  The 
sequences  of  velocity  and  stress  values  at  a  fixed  point  x  in  the 
interior  of  the  target  are  Vq  =  0,  V^; ,  V^,  V|,  Vf , . . .  and 
Z*  =  0.  Z^,  Z*,  Z$,  Z*,..„  with  V\  =  Vx  and  Z^  =  Z,,  as  will  be 
shown  in  the  next  subsection.  On  the  other  hand,  the  sequences 
of  velocity  and  stress  values  at  the  front  face  of  the  target  are 
V\,VL2,VL3,. . .  and  Z,,  Zj,  Z3, . . .;  and  at  the  rear  face  they  are 
Vg  =  0,  V*.  Vr2....  and  Eg  =  0,  Zg,  Eg, . . .. 

The  results  described  above  are  expected  on  physical  grounds.  A 
shock  is  generated  at  the  front  face  of  the  target  at  the  moment  of 
impact:  f  =  t  =  0.  It  travels  at  the  characteristic  (or  acoustic)  wave 
speed  C]  of  the  target  material,  and  arrives  at  the  back  face  of  the 
target  at  t  =  t,  =  l/cu  the  travel  time  (equivalently,  at  the  non- 
dimensional  time  t  =  1 ).  The  shock  reflects  off  the  back  face  and 
arrives  at  the  front  face  of  the  target  at  t  =  2 1,  (equivalently, 
t  =  2).  Thus  the  shock  completes  its  first  round  trip  in  the  target  dur¬ 
ing  the  time  interval  Ij  :  0  sg  t  ^  2t„  (equivalently,  0  ^  r  ^  2).  The 
shock  then  reflects  off  the  front  face  and  completes  a  second  round 
trip  during  the  time  interval  l2  ■  2t,  <  t  s:  4t,  (equivalently, 
2  ^  r  ^  4),  and  so  on  for  all  t  >  0.  The  mth  round  trip  of  the  shock 
wave  (m  =  1 , 2, 3, . . .)  occurs  during  the  time  interval 

Jm  :  2 (m  -  l)t»  C  t  ^  2 mt*.  (98) 


Flyer  t 


Target 


p 

f 


2t, 


t(xL 

N(x,t)  =  0 


Fig.  3.  Regions  of  the  (x,  t)  plane  in  which  the  integer-valued  functions  N  are 
constant.  The  values  indicated  are  assumed  on  the  solid  lines  but  not  on  the  dashed 
lines. 


We  also  will  refer  to  this  interval  itself  as  the  “mth  round  trip”. 
Aside  from  the  instants  2 (m  -  l)t„,  (2m  -  l)t*  and  2 mt,  at  which 
the  shock  is  at  the  boundary  of  the  target,  Im  is  the  union  of  the  open 
subintervals 

C  :2(m-l)t,  <  t  <  (2m  -  l)t„ 

/gj :  (2m  -  l)t*  <  t  <  2mtt. 


8  An  exception  to  the  results  in  this  paragraph  occurs  for  the  special  case  A  =  0. 
Note  that  by  (42)  or  (77),  b  =  0  implies  A  =  0. 

9  This  property  does  not  extend  to  nonlinear  elastic  materials,  for  which  the  shock 
speed  depends  on  the  strength  of  the  shock. 


These  time  intervals  are  illustrated  on  the  right  side  of  Fig.  4.  For 
t  6  (i.e.,  during  the  first  half  of  the  mth  round  trip),  the  shock 

is  propagating  forward  with  position  S^(t)  given  by  (93);  the  region 
of  the  target  to  the  right  of  the  shock  (which  is  ahead  of  it)  has 
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t  t 


Fig.  4.  Locus  of  the  shock  front  in  the  (x,  t)  plane.  M  and  N  are  both  constant  in  the  triangular  shaped  regions,  and  consequently  so  are  the  velocity  and  stress.  VKm  and  £* 
(m  =  0, 1 , 2 _ )  are  the  values  of  the  velocity  and  stress  to  the  right  of  the  shock.  VLm  and  ( rn  1.23 _ )  are  the  values  of  the  velocity  and  stress  to  the  left  of  the  shock. 


velocity  V*  and  stress  E*  ;  and  the  region  to  the  left  of  the  shock 
(which  is  behind  it)  has  velocity  VLm  and  stress  X^.  For  t  e  (i.e., 
during  the  second  half  of  the  mth  round  trip),  the  shock  is  propagat¬ 
ing  backward  with  position  S^(t)  given  by  (93);  the  region  of  the 
target  to  the  right  of  the  shock  (which  is  behind  it)  has  velocity 
V*  and  stress  XjJ,;  and  the  region  to  the  left  of  the  shock  (which  is 
ahead  of  it)  still  has  velocity  VLm  and  stress  Y.Lm,  since  the  reflected 
shock  has  yet  to  alter  this  state.  Finally,  from  the  results  in  the  pre¬ 
vious  subsection  it  follows  that  the  left  values  of  velocity  and  stress 
hold  on  the  forward  propagating  shock,  and  the  right  values  of  veloc¬ 
ity  and  stress  hold  on  the  backward  propagating  shock;  however,  as 
discussed  in  Section  2.3,  values  assigned  to  the  shock  have  no  real 
physical  significance. 

5.3.  Expressions  for  the  velocity  and  stress  on  either  side  of  the  shock 


From  Fig.  4  it  follows  that  the  mth  values  of  the  velocity  and 
stress  on  the  right  side  of  the  shock,  Vjj,  and  X^,,  are  obtained  by 
setting  M(x,t)  =  N(x,t)  =  m  in  the  relations  for  v^(x,  t)  and 
g \  (x,  t)  in  Section  4.6.  In  particular,  on  using  (88)  and  (89)  and 
observing  that  the  impedance  ratios  cancel,  we  obtain  the  follow¬ 
ing  relations  (for  m  =  0, 1,2  . . .): 


VRm  =  V00^-Am),  E*  =  2^(1  -  Am)  =7.2VKm.  (100) 

The  relation  E*  =  z2V*  follows  from  the  fact  that  E^  =  iiV^  (see 
(87)j).  Also,  recall  that  expressions  for  in  terms  of  \f0  or  V i  are 
given  in  (86),  and  expressions  for  E,*,  in  terms  of  V0  or  Xi  are  given 
in  (87).  Similarly,  from  Fig.  4  we  see  that  the  mth  values  of  the 
velocity  and  stress  on  the  left  side  of  the  shock,  VLm  and  E^,  are 
obtained  by  setting  M(x,  t)  =  m  -  1  and  N(x,  t)  =  m  in  the  relations 
for  V\  (x,t)  and  ffi(x,  t)  in  (88)  and  (89).  This  yields  the  following 
relations  (for  m  =  1,2,3,...): 


VL  =  V 

v  m  v 


Z2  ~  Z| 
Zo+Zi 


Am-’ 


M)+Zi 

Vzo  +  z2 


z2  ~Zt 

z2  +z0 


Am 


(101) 


In  deriving  these  relations,  we  used  the  identities 


1 

2 


A  =  z^l 

Zo+Zi 


(103) 


1 

2 


z2  -  Zi  Zp 
Zo+Z,  z2' 


(104) 


The  expressions  on  the  right  in  (101)  and  (102)  follow  from  (86) 
and  (87). 

On  setting  m  =  0  in  (100)  and  using  A0  =  1  (see  (66)),  we  see 
that 


VS  =  E*  =  0,  (105) 

as  indicated  in  Fig.  4.  That  is,  the  region  ahead  of  the  first  forward 
propagating  shock  is  stress-free  and  at  rest,  consistent  with  the 
initial  conditions.  On  setting  m  =  1  in  the  expressions  on  the  right 
(101),  (102)  and  using  A0  =  1,  we  find  that 

V'{  =  V,  and  E^=E,,  (106) 


as  indicated  in  Fig.  4.  That  is,  during  the  first  round  trip  of  the  shock 
the  region  to  the  left  of  the  shock  has  velocity  V i  and  stress  E!  given 
by  (26).  In  particular,  the  conditions  (24)  and  (25)  on  the  front  face 
of  the  target  hold  for  t  <  2 1„  as  claimed  in  Section  3.4. 

On  setting  m  =  1  in  (100),  we  obtain  the  relations  on  the  left 
below  and  also  the  relation  on  the  far  right  for  X^: 


Vo=(l  -A) 


2zi 

Zi  +z2 


V, 


2z0Zi 

(Zo+Zi)(Z,  +z2) 


Vo  >  0, 


(107) 


£?  =  £«,(!  -  A) 


2z2 

Zi  +z2 


2, 


2zoZ]Z2 

(Z0  +  Zi)(Zi  +z2) 


V0  =  z2V«  >  0. 


(108) 


Then  the  expressions  for  and  E^  in  terms  of  V i  follow  from  (86), 
(87)  and  (78)i,  and  the  expressions  in  terms  of  V0  follow  from  (26). 


6.  Additional  properties  of  the  solutions 


E 


L 

m 


(\  Z2  ~  Z1  Zp 

V  z0  +  z,  z2 


=  2, 


Zp  +  Zi  Z2 
Zo  +  z2  Zi 


z2  -  Zi  Zp  ^m_!  \ 
z2  +  z0  z,  )' 
(102) 


Throughout  this  section  we  continue  to  assume  that  z0,Zi,  and 
z2  are  positive  and  finite.  The  limiting  cases  z2  — >  0  or  oo  are 
considered  in  Section  7. 


182 


GA.  Gazonas  et  al.  / International  Journal  of  Solids  and  Structures  75-76  (2015)  1 72-187 


6.1.  Long  time  asymptotes  of  velocity  and  stress 

Let  x  be  any  point  in  the  target.  Recall  that  M(x,  t)  and  N(x,  t)  are 
given  by  (73)  and  (74).  Clearly, 

M(x,  tj— >  oo  and  N(x,  tj^oo  as  t  — >  oo.  (109) 

Since  |A|  <  1,  it  follows  that 

AMm  0  and  Amt)  —  0  as  t  ->  oo.  (110) 

Now  refer  to  the  expressions  (88)  and  (89)  for  the  velocity  and  the 
stress.  On  taking  the  limit  as  t  — >  oo  and  using  (110),  we  see  that 

V\ (x, t)  — >  Voo  >  0  and  0\ (x, t)  — *  I*  >  0  as  t —>  oo.  (Ill) 

Thus  Voo  and  are  the  long  time  asymptotes  of  velocity  and  stress 
in  the  target.  This  conclusion  also  follows  from  (100)-(102).  From 
(86)  and  (87)  we  see  that  Vx  and  2^  depend  only  on  the  impact 
speed  V0  and  the  impedances  of  the  flyer  and  backing  material.  It 
follows  that  at  all  points  in  the  target,  the  long  time  asymptotes  of 
the  velocity  and  stress  are  positive  and  independent  of  the  properties 
of  the  target.10 

Since  the  times  t  in  the  mth  round  trip  Im  of  the  shock  (cf.  (98)) 
approach  oo  as  m  — >  oo,  and  since  jA|  <  1,  from  the  jump  conditions 
in  Appendix  A  we  conclude  that 

[[zq]](t)  — *  0  and  [[or]](t)  — *  0  as  t->  oo.  (112) 

This  result  also  follows  directly  from  (A.4)  and  (110). 

6.2.  Sign  of  the  velocity  and  stress 

During  the  first  round  trip  of  the  shock,  the  velocity  and  stress 
to  the  right  of  the  shock  are  zero  when  the  wave  is  propagating  for¬ 
ward,  that  is  for  0  <  t  <  t,  (see  Fig.  4  or  (105)).  However,  for  t  >  t, 
and  for  all  points  x  to  the  right  of  the  shock,  we  have  V\  (x,  t)  >  0  and 
(x,  t )  >  0,  that  is,  the  particle  velocity  is  in  the  direction  of  impact 
and  the  stress  is  compressive.  Indeed,  for  t  >  t,  the  velocity  and 
stress  to  the  right  of  the  shock  take  on  the  sequence  of  values 
V*  V*,...  and  2*,22, ....  as  determined  by  (100).  Since  |A|  <  1 
and  Voo  and  2^  are  positive,  it  follows  that 

VRm>0  and  2*  >  0  for  m  =  1,2,3, ... .  (113) 

From  (106)  the  velocity  V\  =  Vi  and  stress  2}  =  2i  to  the  left  of  the 
shock  during  the  first  round  trip  are  positive,  a  result  which  is  to  be 
expected  for  the  impact  problem  considered  here.  On  the  other 
hand,  the  sign  of  VLm  and  2}„  for  m  >  1  is  not  so  obvious.  For  those 
cases  where  wave  reflections  would  have  resulted  in  separation  of 
the  flyer  and  target  if  such  separation  had  been  permitted,  we 
expect  that  the“welding”  condition  will  result  in  a  temporary 1 1 
state  of  tensile  stress  at  the  flyer-target  interface  as  well  as  at  all 
points  to  the  left  of  the  shock.  The  problem  that  we  now  address 
is  to  determine  the  conditions  under  which  this  can  occur. 

For  the  remainder  of  this  section  we  assume  that  t  >  2t,  and 
hence  that  t  e  Im  for  some  integer  m  >  2.  The  stress  l,Lm  to  the  left 
of  the  shock  is  given  by  (102),  and  it  follows  that  <  0  if  and  only 
if 

— — —  —  Am~'  >  1.  (tensile  stress  criterion)  (114) 

z0  +  Z\  z2 

Clearly,  this  condition  requires  thatz2  ¥=  Z\  and  A  #  0.  Since  A  can  be 
positive  or  negative  depending  on  the  relative  magnitudes  of  the 
impedances  (see  (77)i),  there  would  seem  to  be  three  cases  to 


consider:  (1 )  z2  >  z,  and  A  >  0  :  (2)z2  >  Zi  and  A  <  0,  in  which  case 
Am~'  must  be  positive,  and  hence  m  odd;  (3)  z2  <  Z\  and  A  <  0,  in 
which  case  Am~’  must  be  negative,  and  hence  m  even.  However, 
we  can  rule  out  the  first  two  cases.  By  (A.3),  2^  =  2®  +  2iAm~1 

for  t  g  /„}.  Since  2j  >  0  and  since  2*_,  >  0  for  m  js  2  (as  shown  in 
the  first  paragraph  above),  we  see  that  a  necessary  (but  not  suffi¬ 
cient)  condition  for  2^  <  0  when  in  ^  2  is  that  A"1-1  <  0,  which  is 
possible  if  and  only  if  A  <  0  and  m  is  even.  This  rules  out  cases  (1) 
and  (2)  above.  Now  restrict  attention  to  case  (3): 

z2<zi,  A  <  0,  me  {2,4,6,...}.  (115) 

We  claim  that  if  2j,  <  0  for  some  n  e  {4, 6 _ },  then  in  fact  I,Lm  <  0 

for  all  m  e  {2 _ ,  n  -  2}.  It  follows  that  if  the  stress  to  the  left  of  the 

shock  ever  becomes  tensile,  then  it  first  becomes  tensile  during  ( the 
entire)  second  round  trip  of  the  shock.  2  To  prove  this,  it  suffices  to 
show  that  if  (114)  holds  for  some  m  =  ne  {4,6,...},  then  it  holds 

for  all  m  e  {2, _ n  -  2}.  But  the  latter  property  follows  easily  from 

the  former  since  |A|m’  >  |A|ra2  for  0  <  mi  <  m2  and  |A|  <1. 

Referring  to  the  expression  (77)i  for  A,  we  see  that  the  inequal¬ 
ities  in  ( 1 1 5 ),  together  with  the  assumption  that  all  impedances  are 
positive,  are  equivalent  to  the  inequalities 

0  <  z2  <  Zi  <  z0.  (116) 


This  and  the  condition  that  m  3=  2  is  even  are  necessary  but  not  suf¬ 
ficient  for  2^  <  0.  A  necessary  and  sufficient  condition  is  obtained 
by  determining  the  additional  restrictions  imposed  on  z0,Zi,z2 
and  m  by  the  inequality  (114).  We  will  restrict  attention  to  the  case 
m  =  2.  From  the  result  in  the  previous  paragraph,  we  know  that  if 
(114)  holds  for  some  even  m  >  2  then  it  is  necessarily  holds  for 
m  =  2.  To  put  this  another  way,  satisfaction  of  (114)  for  m  =  2  is 
necessary  (but  not  sufficient)  for  (1 14)  to  hold  for  other  even  values 
of  m  >  2. 

For  m  =  2,  the  inequality  (114)  reduces  to 


—A  >  1 ; 
Z 0  +  Zi  z2 


(117) 


and  substitution  of  the  expression  (77)j  for  A  yields,  after  some 
algebraic  manipulations,  the  equivalent  inequality 

p(Zio)  <  0,  p(z)  =z2  -  (1  -z20)z  +  3z20,  (118) 


where 


Zi0  =  Zi/Zq,  z20  =  z2/z0.  (119) 

Then  the  inequalities  (116)  are  equivalent  to 

0  <  z20  <  Zio  <  1 .  (120) 

Since  the  graph  of  p(z)  is  concave  up,  p(z)  can  take  on  negative  val¬ 
ues  if  and  only  if  it  has  two  distinct  real  roots,  say  z_  <  z+,  in  which 
case  p(z io)  <  0  if  and  only  if 

z_<z10<z+.  (121) 

The  roots  z±  of  p(z)  depend  on  z20  and  are  given  by 

2z±  =  1  -  2z20  ±  i/q(Z2o),  q(z)=z2-14z  +  l.  (122) 


These  roots  are  real  and  distinct  if  and  only  if  q(z20)  >  0.  The  roots  of 
q(z)  are  7  ±  4^3.  Since  the  graph  of  q(z)  is  concave  up,  q(z20)  >  0  if 
and  only  if  z20  <  7  -  4^3  or  z20  >  7  +  4v/3.  But  the  latter  inequality 
violates  (120),  so  we  must  have 

0  <z20<  7 -4\/3<  0.071797.  (123) 


10  As  will  be  shown  in  Sections  7.1  and  7.2,  this  conclusion  is  also  true  if  the  back 
face  of  the  target  is  either  stress-free  or  rigidly  fixed,  with  the  exception  that  =  0 
in  the  former  case  and  =  0  in  the  latter  case. 

11  That  this  tensile  state  is  necessarily  temporary  follows  from  (111). 


12  If  the  assumption  that  the  flyer  and  target  weld  on  impact  were  not  imposed,  then 
a  tensile  stress  cannot  be  supported  at  the  front  face  of  the  target.  Consequently, 
separation  of  the  flyer  and  target,  if  it  occurs  at  all,  would  necessarily  occur  at 
t  =  tc  =  2t*. 
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Summarizing,  we  have  shown  thatp(Zio)  <  0  if  and  only  if  (123) 
and  (121)  hold,  provided  that  the  middle  and  right  inequalities  in 
(120)  are  also  satisfied  in  this  case.  This  turns  out  to  be  the  case, 
although  we  omit  the  proof.  Then  from  the  earlier  discussion  it  fol¬ 
lows  that  X2  <  0  ( that  is,  the  stress  to  the  left  of  the  shock  is  tensile 
during  the  second  round  trip:  2 1,  <t  <  4 t.)  if  and  only  if  the  impe¬ 
dance  ratio  z20=z2/z0  satisfies  (123)  and  the  impedance  ratio 
z io  =  Zi/z0  satisfies  (121),  with  z±  determined  by  z20  as  in  (122).  Once 
z2 o  and  Zio  have  been  chosen  to  satisfy  these  conditions,  all  three 
impedance  ratios  are  determined  since  Zi/z2  =  zw/z2o:  hence,  the 
selection  ofz20,Zi0  and  one  of  the  impedances  determines  all  three 
impedances.  For  example,  z20  =  l/20  satisfies  (123).  Then 
\fq(z2a)  =  11/20, z_  =  1/5  and  z+  =  3/4,  so  that  (121)  is  satisfied 
if  we  take  zw  =  1/2.  If  we  choose  z,  =  1,  then  z0  =  Zi/Zi0  =  2, 
and  z2  =  z20z0  =  1/10,  and  for  this  set  of  impedances  X2  <  0.  A  plot 
of  the  stress  and  velocity  histories  at  the  impact  face  for  an  impact 
problem  with  this  particular  set  of  impedances  is  given  in 
Section  8. 

7.  Solutions  for  special  cases  of  the  impedance 

In  this  section  we  consider  the  special  cases  z2  =  0  (Section  7.1 ), 
and  z2  =  oo  (Section  7.2).  The  analysis  up  to  this  point  has  also 
assumed  that  z2  is  positive  and  finite,  so  the  results  for  the  cases 
z2  =  0  and  z2  =  oo  are  obtained  by  taking  limits  of  our  earlier 
results  as  z2  approaches  zero  or  infinity. 

7.3.  Target  with  stress-free  back  face  (z2  =  0) 

In  this  section,  we  derive  solutions  for  the  impact  problem 
when  the  back  face  of  the  target  is  stress-free: 

<Ti(/,t)  =  0,  t  e  OS.  (124) 

This  boundary  condition  corresponds  to  the  case  where  there  is  no 
backing  material;  it  also  follows  from  the  back  face  boundary  con¬ 
dition  (32)  in  the  limit  as 

z2  ->  0.  (125) 

Consequently,  the  appropriate  relations  for  the  stress  and  velocity 
in  the  target  for  the  stress-free  boundary  condition  (124)  can  be 
obtained  by  taking  the  limit  of  the  results  in  Sections  4-6  as 
z2  — >  0.  13.  From  (42)  we  see  that  a  is  unchanged  as  z2  — >  0,  whereas 

A^~a  =  zc^  n26) 

Zl  +  Zo 

in  this  case.  By  (85)-(87),  we  see  that  as  z2  — >  0, 

Vl  V0  „  V0  ..  /m, 

1  -  A  ~ *  2  ’  1-4  ’2’  V°°^Vo'  (127) 

On  using  (126)  and  (127)ti2  in  (80)  and  (81),  we  see  that  in  the  limit 
as  z2  — >  0  the  velocity  and  stress  are  given  by 

„l(X)t)=^  [2-4N(v>-4M(w)l, 

2  L  (128) 

<Ti(x,t)=z^[-AN™+AM™], 

for  any  t  >  0  and  0  sg  x  sg  /,  with  A  as  in  (126).  Here  we  have  used 
the  fact  that  the  nonnegative,  piecewise-constant,  integer-valued 
functions  M  and  N  defined  in  (73)  and  (74)  are  independent  of  z2. 


13  Actually,  this  statement  requires  proof,  as  does  the  analogous  statement  for  the 
fixed  back  face  (z2  — *  oo)  in  the  next  subsection.  The  solutions  for  stress-free  or  fixed 
back  faces  can  also  be  derived  directly  by  procedures  analogous  to  those  used  in 
Sections  4-6.  We  have  done  this,  and  the  results  agree  with  those  derived  here  by 
taking  limits  of  the  back  face  impedance.  This  agreement  is  also  corroborated  (for  two 
special  cases)  by  the  discrete  solutions  in  Section  8. 


From  (128)  and  (110),  or  from  (111)  and  (127)34,  we  see  that 

V\  (x,  t)  — >  V0  >  0  and  <7i(x,  t)^0  as  t  — >  oo.  (129) 

At  all  points  in  the  target,  the  long  time  asymptote  of  the  velocity  is 
the  impact  speed  of  the  flyer  and  the  long  time  asymptote  of  the  stress 
is  zero,  independent  of  the  properties  of  the  target. 

Alternative  expressions  for  the  stress  and  velocity  hold  on 
either  side  of  the  shock.  For  m  =  0, 1 , 2, . . .,  the  mth  values  of  veloc¬ 
ity  and  stress  on  the  right  side  of  the  shock  are  obtained  from  (100) 
and  (127)34: 

L*=0,  V*=V0(1  -Am).  (130) 

Note  the  zero  stress  condition  above  is  consistent  with  the  bound¬ 
ary  condition  (124).  For  m  =  1,2,3,...,  the  mth  values  of  velocity 
and  stress  on  the  left  side  of  the  shock  are  given  by 

EIm  =  Z1Am-\  VLm  =  V0(\  --^-Am-^.  (131) 

The  expression  for  VLm  follows  from  (127)3  and  the  relation  on  the 
left  in  (101)  as  z2  — >  0.  The  relation  for  Y,Lm  on  the  left  in  (102)  is 
indeterminate  as  z2  — >  0;  the  relation  on  the  right  approaches  the 
result  above. 

Finally,  by  (131  )i,  we  see  that  the  stress  to  the  left  of  the  shock 
is  tensile  during  the  mth  round  trip  if  and  only  if  A  <  0  and  m  is 
even;  and  by  (126),  A  <  0  iff  z0  >  z,.  Thus  if  z0  >  Zi  then  the  stress 
to  the  left  of  the  shock  is  tensile  during  the  2nd,  4th,  6th,  . .  .round 
trips;  see  Section  8  for  a  particular  example  of  this  case. 


7.2.  Target  with  fixed  back  face  (z2  =  oo) 


In  this  section,  we  derive  solutions  for  the  impact  problem 
when  the  back  face  of  the  target  is  fixed: 

fj(l,t)  =  0.  t  e  M.  (132) 

This  boundary  condition  corresponds  to  the  case  where  the  backing 
medium  is  rigid;  it  follows  from  the  back  face  boundary  condition 
(32)  on  dividing  by  z2  and  then  taking  the  limit  as 

z2  ->  oo.  (133) 


Consequently,  the  appropriate  relations  for  the  stress  and  velocity 
in  the  target  for  the  rigid  boundary  condition  (132)  can  be  obtained 
by  taking  the  limit  of  the  results  in  Sections  4-6  as  z2  — >  oo.  From 
(42)  we  see  that  a  is  unchanged  as  z2  — >  oo,  whereas 

b  — >  1,  A  ^  a  =  (134) 

Zo  +  Zl 

in  this  case.  By  (85)-(87),  we  see  that  as  z2  —  oo, 

Vi  z0  Vo  Zi  „  Vo  n  v  ,1/  / 1  d  c  \ 

1  -  A  ~ *  Zi  2  '  1-A  °2’  (135) 

On  using  (134)  and  ( 1 35)1-2  in  (80)  and  (81),  we  see  that  in  the  limit 
as  z2  — >  oc  the  velocity  and  stress  are  given  by 


Vi(x,t)  =  Z^[-ANM+AM™], 
et, (x,  t)  =  Zo^  [2  -  -  AM(x,t)J . 


(136) 


for  any  t  >  0  and  0  ^  x  ^  l,  with  A  as  in  (134).  From  (136)  and  (110), 
or  from  (111)  and  (135)3f4,  we  see  that 

V\{x,  t)  ^  0  and  03  (x.  t)  — ►  z0V0  >  0  as  t  — >  00.  (137) 

At  all  points  in  the  target,  the  long  time  asymptote  of  the  velocity  is 
zero  and  the  long  time  asymptote  of  the  stress  is  the  product  of  the 
impedance  and  impact  speed  of  the  flyer,  independent  of  the  properties 
of  the  target. 


184 


GA.  Gazonas  et  al.  / International  Journal  of  Solids  and  Structures  75-76  (20 15)  1 72-187 


Alternative  expressions  for  the  stress  and  velocity  hold  on 
either  side  of  the  shock.  For  m  =  0, 1.2 the  mth  values  of  veloc¬ 
ity  and  stress  on  the  right  side  of  the  shock  are  obtained  from  (100) 
and  (135)3f4: 

VRm  =  0,  XRm=z0V0(\-Am).  (138) 

Note  that  the  zero  velocity  condition  above  is  consistent  with  the 
boundary  condition  (132).  For  m  =  1,2,3,...,  the  mth  values  of 
velocity  and  stress  on  the  left  side  of  the  shock  are  given  by 

VLm  =  V,Am-\  ZLm  =  z0V0(\  (139) 

The  expression  for  Y.Lm  follows  from  (135)4  and  the  relation  on  the 
left  in  (102)  in  the  limit  as  z2  — >  oc.  The  relation  for  VLm  on  the  left 
in  (101)  is  indeterminate  as  z2  — »  oo;  the  relation  on  the  right 
approaches  the  result  above. 

Finally,  since  |A|  <  1  and  z0/(z0  +  Zi)  <  1,  from  (138)  and  (139) 
we  see  that  the  stress  is  never  tensile  in  this  case. 

8.  Stress  and  velocity  histories  for  several  impact  problems 

In  this  section  we  plot  the  time  histories  of  the  stress  and 
particle  velocity  at  either  the  front  face  or  the  midpoint  of  the 
target  for  four  different  impact  problems.  The  results  are  shown 
in  Figs.  5-8.  Also  included  in  these  figures  are  the  long-time 
asymptotic  relations  for  the  stress  and  velocity  and,  for  Figs.  5 
and  6,  the  numerical  values  of  these  asymptotes.  In  all  four  exam¬ 
ples  the  velocity  of  the  flyer  is  V0  =  5.  The  Time  on  the  horizontal 
axis  is  the  non-dimensional  time  t  =  t/t„  The  Time  can  also  be 
interpreted  as  the  dimensional  time  t  provided  that  the  travel  time 
t,  =  1,  in  which  case  we  must  have  1  =  Ci  since  t,  =  1/Ci .  The  stress 
and  velocity  values  are  dimensional,  that  is,  04  and  V\. 

Fig.  5  is  a  plot  of  the  stress  and  velocity  histories  at  the  impact 
face  (x  =  0)  for  flyer,  target  and  backing  impedances  of 
z0  =  2,Z]  =  1  and  z2  =  1/10,  respectively.  The  histories  are  calcu¬ 
lated  from  the  relations  (88)  and  (89).  The  impedances  are  taken 
from  the  example  at  the  end  of  Section  6.2,  which  is  a  case  for 
which  the  stress  should  temporarily  becomes  tensile  after  the  first 
round  trip  of  the  shock.  As  is  evident  in  the  figure,  the  stress  is  neg¬ 
ative,  and  hence  tensile,  for  non-dimensional  times  2  <  x  <  4.  Note 
that  since  the  stress  and  velocity  at  the  impact  face  is  not  constant 
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Fig.  5.  Stress  and  particle  velocity  histories  at  the  front  face  of  the  target  (x  =  0)  for 
impact  speed  V0  =  5.  The  flyer,  target  and  backing  impedances  are  z0  =  2,Zi  =1 
and  z2  =  1/10,  respectively.  Observe  that  the  stress  is  tensile  for  non-dimensional 
times  2  <  t  <  4. 
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Fig.  6.  Stress  and  particle  velocity  histories  at  the  midpoint  of  the  target  (x  =  1/2) 
for  impact  speed  V0  =  5.  The  flyer,  target  and  backing  impedances  are  z0  —  4,  z,  =2 
and  z2  =  1.5,  respectively. 
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Fig.  7.  Stress  and  particle  velocity  histories  at  the  midpoint  of  the  target  (x  =  1/2) 
for  a  stress-free  back  face.  The  flyer  and  target  impedances  are  Zo  =  4  and  Z\  =  2, 
and  the  impact  speed  is  V0  =  5. 
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Fig.  8.  Stress  and  particle  velocity  histories  at  the  midpoint  of  the  target  (x  =  1/2) 
for  a  rigid  back  face.  The  flyer  and  target  impedances  are  z0  =  4  and  Z\  =  2,  and  the 
impact  speed  is  V0  =  5. 


I— u  u 

j  9  =  V0zo  T 


- stress  in  target 

- velocity  in  target 

O  discrete  stress  in  target 

*  discrete  velocity  in  target 


»  HI H  I  -*"-.*- 


Ko=0| 


rn 

-e — i  i — e- 


vx  =  Vo  T 


Soc=0| 

-e — »  o  ■ — 


-  stress  in  target 
■  velocity  in  target 
discrete  stress  in  target 
discrete  velocity  in  target 


GA.  Gazonas  et  al.  / International  Journal  of  Solids  and  Structures  75-76  (2015)  172-187 


185 


after  impact,  the  common  “impact-type”  boundary  conditions, 
consisting  of  a  step  in  stress  or  velocity  at  the  front  face  of  the  tar¬ 
get,  would  provide  an  incorrect  solution  to  this  impact  problem. 
Fig.  6  is  for  the  case  of  a  backing  material  with  impedance 
z2  =  1.5;  the  histories  are  calculated  from  the  relations  (88)  and 
(89).  Figs.  6-8  are  plots  of  the  stress  and  velocity  histories  at  the 
midpoint  of  the  target,  x=  1/2,  with  flyer  and  target  impedances 
of  z0  =  4  and  Z\  =  2.  Fig.  7  is  for  the  case  of  a  stress-free  back  face, 
and  Fig.  8  is  for  the  case  of  a  rigid  back  face;  the  time  histories  are 
calculated  from  the  relations  (128)  and  (136),  respectively. 
Since  z0  >  Zi,  the  tensile  stress  during  the  2nd  and  4th  round  trips 
for  the  stress-free  case  in  Fig.  7  is  consistent  with  the  observation 
at  the  end  of  Section  7.1. 

In  a  companion  paper  (Gazonas  and  Velo,  in  preparation),  we 

develop  an  algorithm  for  solving  impact  problems  in  multilayered 
media  with  an  arbitrary  number  of  homogeneous  layers.  The  algo¬ 
rithm  becomes  especially  useful  for  impact  problems  involving 
many  layers,  where  explicit  analytical  solutions  may  be  intract¬ 
able;  see  also  Gazonas  and  Velo  (2012)  and  Velo  et  al.  (2009)  for 
related  work  involving  other  loading  conditions.  The  algorithm  is 
based  on  the  method  of  characteristics  and  involves  a  coupled  sys¬ 
tem  of  recursive  relations  for  the  sequences  of  stress  and  velocity 
values  in  the  interior14  of  each  layer;  a  discrete  version  of  the  impact 
boundary  condition  (18)  is  used.  Closed-form  solutions  of  the  system 
of  recursive  relations  are  obtained  by  using  the  z-transform.  Figs.  6- 
8  include  the  corresponding  discrete  solutions  for  the  particular 
impact  problems  considered  in  these  figures.  The  figures  demon¬ 
strate  the  agreement  between  the  explicit  analytical  solutions  in 
Sections  4.6  and  7  and  the  explicit  discrete  solutions  from  Gazonas 
and  Velo  (in  preparation)  for  all  the  three  back  face  boundary  condi¬ 
tions:  a  semi-infinite  backing  medium,  a  stress-free  or  fixed  back 
face.  In  addition,  for  the  cases  considered  here  we  have  also  verified 
that  the  analytical  and  discrete  solutions  are  in  agreement  with  the 
values  for  the  velocity  and  stress  on  either  side  of  the  shock  obtained 
from  the  expressions  in  Sections  5.3  and  7.15 

9.  Discussion  and  conclusions 

We  begin  with  a  discussion  of  other  methods  (both  classical 
and  more  recent)  for  solving  the  one  dimensional  wave 
equation  in  multilayered,  linear  elastic  media  (Sections  9.1  and 
9.2).  This  is  followed  by  a  discussion  of  the  techniques  used  here 
(Section  9.3). 

9.3.  Time-domain  d'Alembert  methods 

Techniques  based  on  d'Alembert’s  solution  of  the  one¬ 
dimensional  wave  equation  are  quite  common  and  have  been 
extremely  effective;  cf.  Goldsmith  (2001,  p.  38)  or  Graff  (1975,  p. 
95)  for  the  solution  of  simplified  rod  impact  problems  in  linear 
elastic  media.  In  more  recent  studies,  Rossikhin  and  Shitikova 
(2007)  considered  the  thermoelastic  rod  impact  problem  using  a 
time-domain  d’Alembert  formalism.  Similar  time-domain  solutions 
were  developed  by  Bityurin  (2011)  who  investigated  rod  stability 
due  to  repeated  rod  impact  onto  a  rigid  barrier.  Hu  et  al.  (2003) 
considered  the  impact  problem  of  a  rigid  sphere  onto  a  free  rod 
and  validated  a  hybrid  St.  Venant-Hertz  impact  model  with 
experimental  measurements  of  strain  and  velocity  at  fixed 
positions  in  the  rods.  The  problem  of  wave  propagation  induced 


14  For  this  reason,  the  discrete  solution  is  not  plotted  in  Fig.  5. 

15  In  Gazonas  and  Velo  (in  preparation)  we  give  an  algebraic  proof  that  for  the 
impact  problems  considered  in  this  paper  and  for  arbitrary  values  of  the  parameters, 
the  closed-form  discrete  solutions  in  Gazonas  and  Velo  (in  preparation)  are 
equivalent  to  the  the  relations  (100)-(102),  (130)  and  (131)  or  (138)  and  (139)  for 
the  velocity  and  stress  on  either  side  of  the  shock. 


by  impact  was  solved  using  the  d'Alembert  method  by  Yu  et  al. 
(2010)  and  was  subsequently  used  to  verify  an  explicit  space  time 
finite-volume  scheme  for  solving  linear/nonlinear  hyperbolic 
systems. 

9.2.  Other  semi-analytical  methods 

Analytical  methods  other  than  the  d’Alembert  solution  have 
been  used  to  solve  more  difficult  impact  and  wave  propagation 
problems,  and  Laplace  transform  methods  are  almost  universally 
applied  to  transform  and  simplify  the  governing  equations.  In  such 
cases,  a  number  of  either  analytical  or  numerical  inverse  Laplace 
transform  techniques  have  been  used  to  invert  the  transformed 
solutions  back  into  the  time-domain.  For  example,  Schwarz  et  al. 
(2010)  invoked  the  so-called  Laguerre  polynomial  technique  for 
Laplace  transform  inversion  and  derived  solutions  for  rod  interac¬ 
tions  with  deformable  barriers  applicable  to  stamping  tools, 
whereas  Werner  and  Fischer  (1995)  determined  the  inverse 
Laplace  transform  using  the  method-of-residues  for  problems 
involving  the  sudden  arrest  of  rods  in  motion.  The  use  of  residue 
calculus  and  the  evaluation  of  the  Bromwich  integral  (akin  to  the 
Heaviside  expansion  theorem,  cf.  Fodor  (1965,  p.  87))  for  impact 
and  wave  propagation  problems  results  in  solutions  that  exhibit 
the  Gibbs  phenomenon  (Bracewell,  1965),  as  they  are  written  in 
terms  of  an  infinite  series  of  transcendental  functions.  One  may 
also  invoke  the  method  of  separation  of  variables  and  Fourier  anal¬ 
ysis  by  superimposing  waves  of  different  wavelengths  to  generate 
solutions  to  the  wave  equation  (Graff,  1975);  such  methods  of 
superposition  are  generally  inefficient  for  solving  elastodynamic 
problems  with  jump  discontinuities  and  they  also  suffer  from  the 
Gibbs  phenomenon.  More  sophisticated  integral  representation 
solutions  of  the  wave  equation  can  be  derived  using  a  Green’s  func¬ 
tion  formalism  and  often  combine  both  Laplace  and  Fourier  trans¬ 
forms  (Graff,  1975;  Achenbach,  1984). 

Finally,  we  mention  the  class  of  discrete  methods  for  solving  the 
wave  equation  in  multi-layered  media  which,  for  certain  loading 
conditions,  yield  exact  solutions.  These  methods  are  represented 
by  a  set  of  coupled  recursion  relations  for  the  stress  and  velocity 
in  the  layers.  A  variety  of  recursive  solutions  to  the  wave  equation 
have  been  developed  for  impact-type  problems  (Drumheller,  1998; 
Gazonas  and  Velo,  in  preparation),  wave  propagation  and  inverse 
problems  arising  in  geophysics  (Bube  and  Burridge,  1983),  reso¬ 
nance  phenomena  in  multilayered  media  (Gazonas  and  Velo, 
2012),  and  optimal  design  problems  (Velo  et  al.,  2009). 

9.3.  Combined  d'Alembert  and  Laplace  transform  methods 

Less  common  are  solution  methodologies  that  combine  the 
d’Alembert  solution  with  Laplace  transform  methods,  as  illustrated 
in  some  recent  wave  propagation  (Hopkins  and  Gazonas,  2011) 
and  impact  (Randow  and  Gazonas,  2009;  Gazonas  et  al.,  2014; 
Gazonas  et  al.,  2015)  studies.  This  combined  analytical  approach 
was  employed  in  the  current  study. 

Two  additional  features  of  the  current  study  are  the  application 
of  an  impact  boundary  condition  on  the  front  face  of  the  target  and 
the  use  of  the  floor  (or  greatest  integer)  function.  The  impact 
boundary  condition  (1 ),  together  with  the  more  conventional  back 
face  boundary  condition  (32),  allowed  us  to  reduce  the  impact 
problem  to  an  initial-boundary  value  problem  for  the  target  alone. 
Both  of  these  boundary  conditions  are,  with  some  qualifications, 
independent  of  the  target  properties.  The  use  of  the  floor  function 
allowed  us  to  express  the  Heaviside  series  for  the  Laplace  trans¬ 
form  solution  of  the  d’Alembert  functions  in  a  simple  explicit  form. 
This  led  to  apparently  new,  compact  expressions  for  the  stress  and 
particle  velocity  in  the  target  which  are  amenable  to  long-time 
asymptotic  analysis. 


186 


C,.A.  Gazonas  et  al.  / International  Journal  of  Solids  and  Structures  75-76  (20 15)  1 72-187 


Floor  function  constructs  are  used  ubiquitously  in  mathematical 
programming  (Iverson,  1962),  number  theory  (Hardy  and  Wright, 
1979),  and  related  unsolved  number-theoretic  problems  (Guy, 
1994).  Since  a  subset  of  these  solutions  is  exact  (i.e.,  they  exhibit 
no  round  off  error),  they  form  an  important  contribution  to  the 
class  of  benchmark  problems  (Idesman  et  al.,  2009)  relevant  to  ver¬ 
ification  of  large-scale  computational  codes  (Oberkampf  and  Roy, 
2010). 

Our  results  are  corroborated  by  discrete  solutions  of  a  recursive 
system  of  equations,  derived  in  Gazonas  and  Velo  (in  preparation), 

and  applied  to  the  impact  boundary  value  problem  studied  here. 
On  the  other  hand,  the  agreement  of  the  two  methods  provides 
partial  verification  of  the  these  recursive  methods  for  this  special 
case.  This  is  important  since  exact,  analytical  solutions  for  a  wider 
class  of  impact  problems  involving  targets  with  multiple  layers  are 
difficult  to  obtain. 

Our  current  work  in  this  area  indicates  that  the  impact  bound¬ 
ary  condition  (1)  is  also  applicable  to  the  impact  of  a  semi-infinite 
elastic  flyer  on  a  piezoelectric  target  backed  by  a  semi-infinite  elas¬ 
tic  half-space  (Gazonas  et  al.,  2014).  Furthermore,  this  boundary 
condition  can  be  modified  to  apply  to  both  elastic  and  piezoelectric 
flyers  of  finite  thickness  (Gazonas  et  al.,  2015). 

The  use  of  floor  functions  for  solving  problems  in  elastodynam- 
ics  has  implications  beyond  problems  involving  impact.  For  exam¬ 
ple,  the  classic  problem  of  a  finite  elastic  strip  fixed  at  x  =  0  with  a 
transient  traction  p(t)  applied  at  x  —  I,  cf.  Eringen  and  §uhubi 
(1975,  Section  6.6),  can  also  be  written  in  explicit  form  using  floor 
functions.  Also,  the  solution  for  a  string  (with  unit  wave  speed) 
subjected  to  a  transient  displacement  f(t)  at  x  =  0  with  the  end 
x  =  I  allowed  to  slide  freely  along  a  line  perpendicular  to  the  unde¬ 
formed  string,  cf.  Courant  and  Hilbert  (1962,  pp.  508-510),  is  also 
solvable  in  terms  of  floor  functions. 


Appendix  A.  The  jumps  in  stress  and  velocity  across  the  shock 

Let  S(t)  denote  the  position  of  the  shock  at  time  t  ^  0.  At  an 
instant  t  when  S(t )  is  in  the  interior  of  the  target,  that  is, 
0  <  x  <  I,  the  jump  in  a  field  variable  tp  across  the  shock,  denoted 
by  [[(p]],  is  defined  as  follows: 

M\{t)  = 

frombehindtheshock  fromaheadoftheshock 

=  lim tp(S(t),  t)  -  lim cp(S(t),t).  (A.l) 

tu  Ht 

The  opposite  sign  convention,  that  is,  (A.l)]  with  the  two  limits 
interchanged,  is  also  used  in  some  of  the  literature.  The  equivalence 
of  these  relations  is  easily  seen  from  Fig.  4  if  we  recall  that  for  a  for¬ 
ward  propagating  shock,  positions  x  to  the  right  of  the  shock  are 
ahead  of  it  and  positions  x  to  the  left  of  the  shock  are  behind  it, 
while  for  a  backward  propagating  shock  the  reverse  is  true.  The 
relation  on  the  right  is  also  well-defined  at  an  instant  t  when  S(t) 
is  on  the  front  or  back  face  of  the  target,  so  this  is  the  definition 
of  the  jump  that  will  be  used  in  these  cases.  Observe  that  [[<p]]  is 
positive  (negative)  if  and  only  if  the  value  of  cp  behind  the  shock 
is  larger  (smaller)  than  the  value  of  tp  ahead  of  the  shock,  equiva¬ 
lently,  if  tp  increases  (decreases)  on  the  passage  of  the  shock.  Of 
course,  [[<p]]  =  0  if  and  only  if  cp  is  continuous.  The  shock  wave  is 
compressive  (tensile)  at  the  instant  t  if  [[<7i]](t)  is  positive  (negative), 
that  is,  if  the  stress  behind  the  shock  exceeds  (is  exceeded  by)  the 
stress  ahead  of  it.  Note  that  this  definition  refers  to  the  sign 
of  the  jump  in  stress,  not  to  the  signs  of  the  stress  on  either  side 
of  the  shock.  Now  consider  the  jump  in  velocity  and  stress  across 
the  shock  during  the  mth  round  trip  Im  (m  =  1, 2, 3, . . .).  Restricting 
attention  to  those  instants  t  at  which  the  shock  front  S(t)  lies  in  the 
interior  of  the  specimen,  we  claim  that 
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The  relations  on  the  left  follow  immediately  from  Fig.  4  and  the  def¬ 
inition  (A.l)  of  the  jump.  The  top  expressions  on  the  right  can  be 
obtained  from  the  relations  (TOO)— (102)  for  V^.Z*.  VLm  and  HLm  after 
some  algebraic  manipulations;  we  can  also  proceed  as  in  the  next 
paragraph.  The  bottom  expressions  on  the  right  follow  from  (77)2 
and  (26).  From  (80)  and  (81)  we  see  that  the  jumps  in  velocity 
and  stress  at  any  time  t  ^  0  are  given  by 

[[^i]](t)=/Lr[[[^]](t)-fa[[>lM]](t)], 

71  “  1  L  J  (A.4) 

[MKO^I^Ktj  +  b^Kt)]. 

Consider  the  case  where  t  e  /„(.  Then  from  Figs.  2,  3  or  4,  M(x,  t)  is  a 
continuous  function  of  x  and  t,  and  (A.4)  reduces  to 

[MK0  =  ^  [[ANm,  Mt)  =  [[AN]](t). 

From  Fig.  4  and  the  definition  (A.l)  of  the  jump,  we  see  that 
[[Aw]](t)  =  Am  -Am~’  =  (A-  l)Am~\  which  together  with  the  above 
relations  yields  the  top  middle  expressions  in  (A.2)  and  (A.3).  Sim¬ 
ilar  arguments  yield  the  bottom  middle  expressions;  note  that  in 
this  latter  case  N(x,  t)  is  continuous  and  M(x,  t)  suffers  a  jump  dis¬ 
continuity.  This  completes  the  proof  of  (A.2)  and  (A.3). 

From  (42)  we  see  that  the  sign  of  A  depends  on  the  signs  of  the 
impedance  differences  z0  -  Z\  and  z2  —Z\.  Consequently,  by  (A.2) 
and  (A.3),  the  signs  of  the  jumps  in  velocity  and  stress  depend  on 
the  signs  of  these  impedance  differences  and,  if  A  <  0,  on  whether 
m  is  even  or  odd.  From  the  middle  or  left  expressions  in  (A.2)  and 
(A.3),  we  see  that  the  jump  in  velocity  (stress)  across  the  backward 
propagating  shocks  is  a  constant  multiple,  namely  - b  (+b),  of  the 
corresponding  jump  across  the  forward  propagating  shocks.  Also, 
from  the  middle  expressions  for  [[iq]](t)  and  [[or]](t)  and  the  fact 
that  Z]  =Z]V]  (see  (26)),  or  directly  form  the  expressions  on  the 
right  in  (A.2)  and  (A.3),  we  see  that 

[[(Ti]](t)  =  ±Zi[[zq]](t)  (in  the  interior  of  the  target),  (A.5) 

independent  of  m  and  also  independent  of  the  properties  of  the  flyer 
and  backing  medium.  Here  the  “+”  case  holds  for  forward  propagat¬ 
ing  shocks  and  the  ”  case  holds  for  the  backward  propagating 
shocks.16 

From  (66)  and  (A.2)-(A.3)  with  m  =  1,  or  directly  from  Fig.  4,  we 
see  that  the  jumps  in  velocity  and  stress  across  the  first  forward 
propagating  shock  (t  e  Ifi  :  0  <  t  <  t ,)  are  given  by  V ]  and  Zj, 
respectively.  These  are  also  the  jumps  in  velocity  and  stress  at 
the  front  face  of  the  target  at  the  moment  of  impact,  t  =  0.  How- 


16  The  relation  (A.5)  can  be  established  independently  of  (A.2)  and  (A.3).  Since  there 
is  no  wave  reflection  during  either  of  the  time  intervals  in  question,  momentum 
balance  across  the  shock  is  given  by  [[<Ji]]  =  f>]  ii,  [[iq]],  where  Ut  is  the  shock  velocity 
(i.e.,  the  signed  speed):  LI,  =  +c, ,  with  the  +  (-)  case  holding  when  the  wave  in  the 
target  is  propagating  in  the  positive  (negative)  x  direction. 
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ever,  the  definition  of  the  jump  on  the  left  side  of  (A.1)  does  not 
apply  at  the  front  or  back  face  of  the  target  for  times  t  >  0,  and 
the  relations  (A.2)  and  (A.3)  are  not  valid  there.  At  the  front  and 
back  faces,  the  velocity  and  stress  suffer  jumps  at  any  instant  when 
shock  reflection  occurs,  and  we  apply  the  expression  for  the  jump 
on  the  right  in  (A.l).  For  the  back  face,  shock  reflection  occurs  at 

t  =  t,,  3t„,  5t . (see  Fig.  4),  and  the  jumps  in  velocity  and  stress 

at  these  instants  are  given  by 

[[*'i]]((2m  -  l)t.)  =  VRm  -  =  V*Am  \  (A.6) 

M(2m  -  l)t.)  =  Z*  -  ZRm_,  =  (A. 7) 

for  m  =  1,2,3,....  Recall  that  expressions  for  Vf  and  2,  are  given  in 
(107)  and  (108).  For  the  front  face,  shock  reflection  occurs  at 
t  =  2t„,4t,,6t,, . . .,  and  the  jumps  in  velocity  and  stress  at  these 
instants  are  given  by 

[MP mt.)  =  Ci  -  <  =  •  Cm-' ,  (A.8) 

z  0  +Z\ 

[M(2l7lt.)  =  Cl  -  2m  =  rTTT'  ^ra-1,  (A.9) 

for  m  =  1,2,3, . . ..  Note  that  these  are  constant  (i.e.,  independent  of 
m)  multiples  of  the  corresponding  jumps  on  the  back  face.  The  rela¬ 
tions  on  the  left  in  (A.6)-(A.9)  follow  easily  from  Fig.  4.  The  relations 
on  the  right  follow  from  (100)— ( 1 02)  and  ( 1 07)— ( 1 08).  They  can  also 
be  obtained  from  (A.4);  note  that  M  and  N  suffer  jump  discontinu¬ 
ities  simultaneously  in  these  cases.  Since  ZR  =  z2VR  (see  (108)),  the 
relations  (A.6)-(A.9)  imply1'  that  for  t  =  t,,3t,,  5t„ . . ., 

[[<7i]](t)  =  z2[[z'i]](t)  (on  the  back  face);  (A.10) 

and  for  t  =  2t,,4 1,,  6 1„, . . ., 

[[cTi]](t)=-Zo[[zh]](t)  (on  the  front  face).  (A.l  1) 
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